m 

CM 

00 

< 

I 

a 

< 


MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
RESEARCH  LABORATORY  OF  ELECTRONICS 
CAMBRIDGE,  MASSACHUSETTS  02139 


Final  Report 
covering  the  period 
1  August  1981  -  31  May  1984 

Sponsored  by 

Advanced  Research  Projects  Agency 
ARPA  Order  No.  4258 


Program  Code  Number: 

80230 

Contractor: 

Massachusetts  Institute  of  Technology 


Contract  Number: 

N0001 4-81 -K-0662 . 

/ 

Principal  Investigator: 
Jeffrey  H.  Shapiro 


Scientific  Officer: 

Director  Earth  and  Environmental  Physics  Program 

Arctic  and  Earth  Sciences  Oi vis  ion 

Office  of  Naval  Research 

800  North  Quincy  Street 

Arlington,  Virginia  22217 

Title  of  Work: 


Two-Photon  Coherent  State  Light  —  Its  Generation 
and  Potential  Applications 


'^!aEC  6  108*^^ 


Effective  Date  of  Contract: 
August  1 ,  1981 


Contract  Expiration  Date: 
May  31,  1984 


Amount  of  Contract: 
$300,000 


„  nlrtVplC  ^  • 


itsVA- 


riu  CORY 


•  '  *  ' 


w\: 

Jt 


•** 1  ;  r  I  ' 


1 


26  066 


84  n 


Two-Photon  Coherent  State  Light  — 

Its  Generation  and  Potential  Applications 


13 


Abstract  13 

Research  Summary  14 

References  24 

Personnel  27 

Publications  28 

Appendix  I 
Appendix  II 
Appendix  III 


'-'in  n. ... 

•  v/l\  K  '■ 


•  f-  RELEASE 
’‘<•’1  TED 


A^onsion  Foi; 


□afr 


-  13  - 
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Abstract 

\ 

Recent  work  has  highlighted  the  potential  applications  of  two- 
photon  coherent  states  (TCS),  also  known  as  squeezed  states,  in  optical 
communi cations  and  precision  measurements.  This  research  program  was 
aimed  at  generating  TCS  light  via  degenerate  four-wave  mixing  (DFWM), 
and  verifying  the  non-classical  nature  of  TCS  light  via  photon¬ 
counting  measurements.  The  preceding  experimental  work  was  supported 
by  analytical  studies  of  DFWM  TCS  generation,  and  quantum  photodetection. 
In  addition,  analyses  were  performed  of  key  proposed  TCS  applications  in 
optical  communi cations  and  phase-sensing  interferometry.  Although  the 
experimental  effort  did  not  yield  any  TCS  observations,  it  did  succeed 
in  providing  the  first  quantum-noise  limited  measurements  obtained  from 
DFWM.  Moreover,  the  associated  theoretical  effort  clarified  pump  quantum 
noise,  probe-conjugate  loss,  and  backward  vs  forward  DFWM  issues  in  TCS 
generation.  Furthermore,  the  applications  research  developed  simultane¬ 
ous  amplitude  and  phase  uncertainty  relations  for  optical  heterodyne 
detection,  and  showed  explicitly  how  to  use  TCS  to  surpass  the  so-called 
standard  quantum  limit  (SQL)  on  position  sensing  in  a  gravity-wave 
detecting  interferometer. 


RESEARCH  SUMMARY 


A  highly  stabilized  laser  produces  a  light  beam  in  a  Glauber  coherent 
state.  This  means  the  beam  is  endowed  with  the  following  quantum  fluctua¬ 
tion  behavior:  photon  counting  measurements  are  Poisson  distributed  hence 
the  count  variance  equals  the  count  mean;  and  homodyne  measurements  are 
Gaussian  distributed  with  signal-to-noise  ratio  equal  to  four  times  the 
average  number  of  received  photons.  Furthermore,  the  observed  fluctuation 
behavior  for  any  laser,  light-emitting  diode,  or  ordinary  light  source  can 
be  accounted  for  by  averaging  the  preceding  results  over  a  classical  ensemble 
representing  the  light  beam  as  a  random  mixture  of  coherent  states.  These 
results  form  the  foundation  for  all  of  the  semi  classical  optical  detection 
and  communication  analyses. 

Recent  theory  has  shown  that  the  generation  of  light  beams  with  quantum 

states  of  superior  fluctuation  behavior,  relative  to  Glauber  coherent  states, 

may  be  possible.  These  states  are  called  two-photon  coherent  states  (TCS), 

or  squeezed  states.  They  are  basically  minimum  uncertainty  states  for  the 

quadrature  components  of  the  electromagnetic  field  possessing  an  asymmetric 

noise  division  between  the  quadratures,  so  that  one  quadrature  has  a  "squeezed" 

noise  variance,  i.e.,  it  is  less  than  the  coherent  state  value  of  1/4.  In 

photon  counting  measurements  at  high  mean  counts  TCS  light  can  have  a  count 

2/3 

variance  as  low  as  N$  '  ,  where  Ns  is  the  average  number  of  received  photons. 

In  homodyne  measurements,  because  of  the  quadrature  noise  squeezing,  TCS 
light  can  have  a  signal-to-noise  ratio  as  high  as  4NS(N$+1). 

Under  this  research  program,  we  have  undertaken  an  experimental  effort 
aimed  at  generating  TCS  light,  and  verifying  some  of  its  novel  fluctuation 
characteristics  [l]-[5].  In  concert  with  the  experimental  activity,  we 
have  worked  toward  improving  the  theory  behind  TCS  generation  schemes  [5] -[7], 


advancing  understanding  of  quantum  photodetection  (especially  as  applied  to 
TCS  light)  [8],.  [9],  and  analyzing  potential  applications  of  TCS  in  optical 
communications  [10]  and  phase-sensing  interferometry  [11],  [12]. 

TCS  Generation  -  Theory 

Two-photon  coherent  states  may  be  generated,  in  principle,  by  a  variety 
of  nonlinear  optical  processes,  including  the  degenerate  parametric  amplfier 
(DPA)  and  degenerate  four-wave  mixing  (DFWM).  The  essence  of  generating 
TCS  may  be  stated  simply  as  follows.  For  a  single-mode  field  with  annihila- 
tion  operator  a,  mix  a  part  of  this  field  with  its  phase  conjugate  field 
represented  by  the  creation  operator  a  to  create  a  new  mode  represented  by 

jm. 

the  operator  c  such  that 

c  =  ua  +  vaf  ,  ( II. 1 ) 

where  |y|  -  |v|  =1  ensures  that  c  is  an  annihilation  operator.  Then,  if 

A  A 

the  mode  a  is  in  a  coherent  state  (CS),  the  mode  c  will  be  in  a  TCS.  Thus, 
a  physical  process  that  generates  a  phase  conjugate  field  for  some  input 
field  is  a  possible  candidate  for  generating  TCS  light.  DFWM  is  such  a 
process,  and  was  suggested  by  Yuen  and  Shapiro  [13]  as  a  possible  source  of 
TCS.  In  the  standard  backward  DFWM  geometry  treated  in  [13],  two  counter- 
propagating  pump  waves  intersect  an  object  (probe)  wave  at  a  small  angle  in 
a  nonlinear  medium  possessing  a  third  order  (x^)  nonlinearity.  All  three 
waves  are  at  the  same  frequency.  A  phase-matched  interaction  in  the  non¬ 
linear  medium  generates  a  phase-conjugate  (PC)  image  wave  which  propagates 
in  the  opposite  direction  to  the  object  wave.  The  outputs  of  the  four-wave 
mixer  can  then  be  separated  by  means  of  isolators,  and  combined  through  an 
optical  delay  line  with  the  proper  phase  relationships  on  a  50/50  beam 


splitter.  In  a  classical  analysis,  the  complex  field  amplitudes,  and  A^, 


of  the  image  and  object  waves  obey  the  coupled-mode  equations 


dA, 


*  ★ 


dz  1K  A2  » 


(II. 2a) 


dA. 


dz 


★  * 

•i< 


(II. 2b) 


*  -1  (3) 

where  k  =  2  u>x  B-^/cn  gives  the  complex  coupling  constant  in  terms  of  the 
nondepleted  pump  waves  of  complex  field  amplitudes  B.|  and  82*  Yuen  and 
Shapiro  [13]  gave  a  quantum  version  of  the  above  treatment  assuming  that  the 


pump  fields  B-|  and  B2  are  strong  and  can  remain  classical  along  with  the 


(3) 


medium  described  by  the  third-order  nonlinear  susceptibility  x  •  They  quan¬ 


tized  the  object  and  the  image  modes  replacing  A.  and  A.  with  the  photon 

J  J 


annihilation  and  creation  operators  a.  and  a-+,  respectively,  for  j  *  1 ,  2. 

J  J 


Equation  (I I. 2)  then  becomes 


da. 


★« 


dz 


*  Ik  a. 


(II. 3a) 


^a2  *-  + 
d r  3  ai 


(II. 3b) 


with  the  following  solutions 

a1(0)  =  ua-| (L)  -  iva2+(0)  ,  (II. 4a) 

a2(L)  »  ya2(0)  -  1va  +  (L)  ,  (II. 4b) 


where  y  =  sec(|<|L),  v  3  e"1 9tan(  |k |L) ,  <  =  |<|ei0,  a^L)  and  a2(0)  are  the 
input  fields  to  the  four-wave  mixer  at  z  =  L  and  z  =  0,  respectively. 

When  the  outputs  of  the  four-wave  mixer  are  combined  through  a  50/50 
beam  splitter  to  generate  two  new  modes 


c  =  [3,(0)  -  ia2(l)]/21/2  ,  (II. 5a) 

d  =  [a,(0)  +  ia,(L)]/21/2  ,  (II. 5b) 


the  solutions  become 


uc 


in 


-  -  -  X 

d  *  ud.  -  vd. 
in  in 


in  terms  of 

c1n  3  [a-j (L)  -  ia2(0)]/21/2  , 
din  3  [a-j ( L)  +  ia2(0)]/21/2  , 


(II. 6a) 
(II. 6b) 

(II. 7a) 
(II. 7b) 


which  are  linear  combinations  of  the  input  modes  to  the  four-wave  mixer. 

2  2 

Equation  (II. 6)  is  of  the  same  form  as  (II. 1)  and  j u |  -  |v[  =  i ,  so  that 

A  A  A  A 

the  modes  c  and  d  are  in  TCS  if  a^(L)  and  a2(0)  are  in  CS. 

The  above  analysis  indicates  that  backward  DFWM  is  a  source  of  pure  TCS. 
Depending  upon  the  phase  and  magnitude  of  v,  arbitrary  noise  squeezing  is 

A  A 

predicted  to  be  obtained  in  one  of  the  quadratures  of  mode  c  or  mode  d.  In 
a  real  experiment,  this  is  not  quite  so  because  of  the  assumptions  made  in 
arriving  at  (II. 3).  In  particular,  the  pump  modes  B-j  and  B2  cannot  neces¬ 
sarily  be  considered  classical,  and  the  effects  of  their  quantum  amplitude 

A  A 

and  phase  fluctuations  on  modes  c  and  d  should  be  calculated.  Moreover,  the 
preceding  analysis  neglects  loss,  which  is  known  to  have  a  strong  effect  on 
quantum  statistics. 

In  [5],  [6]  we  have  reported  our  analyses  of  generalizing  the  Yuen  and 
Shapiro  model  to  encompass  pump  quantum  noise  in  a  lossless  backward  DFWM 
setup,  and  loss  on  the  probe  and  conjugate  beams  in  backward  DFWM  in  the 
absence  of  pump  quantum  noise.  Both  of  these  studies  focused  on  the  quadra- 

A  A 

ture  noise  squeezing  in  the  output  modes  c  and  d  from  (1 1. 6).  In  particular, 
for  the  quadrature  variance  <Ac^>  =  <A  [(c+c  )/2\c>  it  was  shown  that  the 
Yuen  and  Shapiro  result  is  obtained  in  a  lossless  case  with  pump  quantum 
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noise  included  if  the  pump  beams  are  very  intense  ( |B^  |  =  jB^J  — ®)  the 
coupling  is  very  weak  (2’^wx^/cn  — *■  0)  and  the  gain  is  constant 
( | k | L  *  constant).  Thus,  pump  quantum  noise  places  no  fundamental  limitation 
on  TCS  generation. 

-  2 

The  results  obtained  for<Ac-|  >  in  backward  DFWM  with  probe-conjugate 
loss  were  not  so  favorable.  Here  it  was  found  that  even  a  moderate  amount 
of  such  loss  (e.g.,  loss  coefficient  y  =  0.5/L)  severely  limited  the  obtain¬ 
able  noise  squeezing.  Thus,  loss  presents  a  severe  restriction  on  the  use  of 
resonant-medium  backward  DFWM  as  a  TCS  generator. 

In  an  attempt  to  circumvent  the  preceding  limitation,  we  turned  our 
attention  from  backward  DFWM  to  forward  DFWM  [7].  In  forward  DFWM,  two 

(31 

strong  pump  beams  propagate  at  small  angles  ±  $/2  from  the  +z  axis  in  a  x 
medium.  A  phase  matched  interaction  ensues  between  probe  and  conjugate  waves 
at  angles  ±  <j>/2  from  the  +z  axis  in  a  plane  obtained  by  rotating  the  pump- 
beam  plane  through  an  arbitrary  angle  about  the  z  axis.  This  nominally  co¬ 
propagating  geometry,  when  analyzed  quantum  mechanically  along  the  lines  of 
the  Yuen  and  Shapiro  backward  DFWM  case,  yields 

a-j  (L)  =  ya-j  (0)  -  iva2+  (0)  ,  (II. 8a) 

a2+(L)  =  ya2+(0)  +  iv*a.,(0)  ,  (II. 8b) 

for  the  output  annihilation  operators  at  z  =  L  in  terms  of  the  input  annihila- 

tion  operators  at  z  *  0,  where  y  =  cosh  (|<|L),  v  =  e’  sinh  ()<]L),  and 
i  9 

<  =  |<|e  .  Defining  new  output  and  input  modes  via  (cf.  (II. 6),  (II. 7)) 

c  =  [a1 (L)  -  ia2(L)]/21/2  ,  (II. 9a) 

1/2 


d  =  [a-j(L)  +  ia2(L)]/2 


(II. 9b) 
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and 

cin  s  [ai (0)  -  ia2(0)]/21/2  , 

dfn  a  [a^O)  +  ia2(0)]/21/2  , 
yi el ds 

A  *  A  I, 

c  =  uC ■  -  VC. 

*  in  in  * 

and 

A  A  A  «, 

d  =  ud-  +  vd-  . 

nn  in 

Thus,  for  lossless  forward  DFWM  with  no  pump  quantum  noise  ideal  TCS  genera¬ 
tion  results  from  50/50  interferometric  combination  of  the  two  output  beams. 

We  have  not  explicitly  examined  the  pump  quantum  noise  behavior  of  for¬ 
ward  DFWM;  we  expect  it  will  mimic  that  of  backward  DFWM.  We  have  shown, 
however,  that  probe-conjugate  loss  does  not  present  a  fundamental  limit  on 
obtainable  quadrature  noise  squeezing  from  forward  DFWM.  This  is  because 
the  counter-propagating  geometry  of  backward  DFWM  causes  a  looping  of  loss- 
induced  quantum  noise  in  the  interaction  medium  that  is  absent  in  the 
co-propagating  geometry  of  forward  DFWM.  Thus,  forward  DFWM  is  a  more  promis 
ing  TCS  generator  than  is  backward  DFWM. 

TCS  Generation  -  Experiment 

We  chose  to  pursue  TCS  generation  experimentally  via  backward  DFWM  in 
sodium  vapor.  (This  work  was  initiated  before  our  loss  analysis  was 

(31 

completed.)  Because  DFWM  occurs  through  the  third-order  nonlinearity  x 
of  the  mixing  medium,  we  employed  a  pulsed  optical  system  and  exploited  the 
resonant  enhancement  of  sodium  xv  in  the  vicinity  of  539  nm  wavelength  D2 
transition  in  order  to  access  the  high  reflectivity  regime.  Furthermore, 
because  a  pulsed  experiment  did  not  lend  itself  to  homodyne  detection,  we 


(II. 10a) 
(II. 10b) 

(II. 11a) 

(II. lib) 


used  photon  counting  to  look  for  anti-bunching  (i.e.,  second  reduced  factorial 


moment  g£  less  than  unity)  in  light  that  had  undergone  DFWM.  No  previous  DFWM 
experiments  had  ever  tried  to  look  at  quantum  noise;  it  turned  out  that  our 
initial  DFWM  quantum-noise  experiment  was  overconstrained  by  a  combination  of 
factors,  as  elaborated  below. 

The  experiment  comprised  three  subsystems:  a  continuous -wave  oscillator/ 
pulsed-amplifier  chain  dye  laser,  a  heat-pipe  oven  sodium  vapor  cell,  and  a 
special  purpose  photomultiplier  tube/computer  data  acquisition  photon  counting 
setup.  The  laser  system  [2]  was  developed  to  approach  insofar  as  possible  the 
stable  center- frequency  transform-limited  pulse  train  needed  for  eliminating 
classical  excess  fluctuations  in  the  DFWM  conjugate  output  beams  [4].  The 
photon  counting  system  [1]  was  assembled  to  permit  accurate  photon  statistics 
determination  from  the  nanosecond  duration  pulses  produced  by  the  laser.  As 
described  in  [3],  [5]  background  light  from  the  sodium  cell  had  to  be  reduced 
to  less  than  one  detected  photon  per  pulse,  by  means  of  spatial  filtering, 
time  gating,  and  adjustment  of  DFWM  pump  strength  and  pump/probe  angular  sepa¬ 
ration,  in  order  to  make  the  desired  quantum  noise  measurements.  The  result 
was  that  we  had  to  operate  in  a  low- reflectivity  regime  for  which  Poisson 
statistics  (g2=l)  were  expected.  Our  measurements  confirmed  the  Poisson 
behavior  [3],  [5],  indicating  that  we  were  able  to  reduce  excess  fluctuations 
to  acceptably  low  levels,  thus  demonstrating  the  first  quantum- limited  noise 
measurements  obtained  on  a  DFWM  output  beam. 

Quantum  Photodetection 

As  an  adjunct  to  the  detection  aspects  of  the  TCS  generation  work,  and 
as  a  foundation  for  TCS  applications  research  in  optical  communications  and 
precision  measurements,  a  number  of  new  theoretical  results  were  developed 
in  the  area  of  quantum  photodetection.  Building  from  Cook's  photon  flux 
picture  of  photodetection  [14],  we  relaxed  [8],  [9]  the  quasimonochromati c 


condition  in  Yuen  and  Shapiro's  operator  representations  [15]  for  direct 
detection, .homodyne  detection,  and  heterodyne  detection,  and  extended  [8]  the 
coherent  detection  results  of  [15]  to  include  dual  detection  configurations 
(cf.  [16])  and  local  oscillator  excess  noise.  Moreover,  the  results  of  [8] 
were  .couched  in  terms  which  made  comparison  with  familiar  multi -mode  semi- 
classical  shot  noise  particularly  simple.  Finally,  in  [9],  we  showed  that 
multi -mode  TCS  permit  heterodyne  detection  amplitude  and  phase  measurements 
to  be  made  simultaneously  without  being  subject  to  any  uncertainty  principle. 
This  phenomenon  has  applications  in  precision  measurements  (see  below). 

Optical  Communications  with  TCS 

Helstrom  [17],  [18]  found  the  optimum  quantum  measurement  operator  for 
deciding  between  equally  likely  density  operators  and  for  the  state  of 
a  quantized  signal  field  to  be  as  follows.  The  minimum  error  probability 
decision  rule  chooses  hypothesis  H.  ( i . e . ,  state  p.)  when  the  outcome  of  mea- 

J  J 

suring  u(p-j  -  pQ)  is  j,  where  u(x)  is  the  unit  step  function.  Helstrom  also 
showed  that  the  error  probability  of  this  receiver  obeys 

Pr(e)  =  Z‘\l  -  [1  -  |<,lIv|2]1/2>  .  (H-12) 


for  the  pure  state  problem  p  .  =  |«f  .><'?.•.  |.  Kennedy  [19]  and  Dolinar  [20]  later 

J  J  J 

deduced  realizations  for  exponentially  optimum  and  exactly  optimum  receivers, 
respectively,  when  the  \v  ■>  are  coherent  states,  using  conditional  Poisson 
process  photodetection  models.  Shapiro  [21],  using  the  representation  theorem 
from  [15],  was  able  to  explicitly  develop  the  operator  description  of  the  near¬ 
optimum  receiver.  In  this  program  we  have  extended  the  work  of  Kennedy,  Dolinar, 
and  Shapiro  to  the  detection  of  binary  two-photon  coherent  state  signals. 


In  comparison  with  coherent  state  signals  TCS  offer  significant  communi¬ 


cation  performance  advantages,  e.g.,  (11.12)  reduces  to 


r 


Pr(e)  =< 


2''{1  -  [1  -  exp(-4Ns)]1/2}  , 

for  coherent  state  |*.> 

J 


(11.13) 


2"1{1  -  [1  -  exp(-4N$(Ns+l))]1/2}  , 

^  for  optimized  TCS  |h'j>  , 


of  average  energy  content  Ns  photons.  We  have  developed  a  general  explicit 
realization  for  the  u(p^  -  pg)  measurement  for  TCS  signals  [10].  Our  construc¬ 
tion  of  the  optimum  TCS  measurement  combines  the  operator  transformation  prop¬ 
erties  of  the  idealized  four-wave  mixer  in  [13]  first  with  the  near-optimum 
operator  analysis  in  [21]  and  later  with  a  quantum  version  of  the  Dolinar 
work  [20]. 


Phase-Sensing  Interferometry 

The  performance  of  phase-sensing  interferometers  employing  TCS  and  homo¬ 
dyne  detection  were  analyzed  [11],  [12]  and  compared  to  the  performance  of 
systems  employing  direct  detection  [22].  Standard  differenced  direct-detection 
Michel  son  and  Mach-Zehnder  interferometers  were  shown  to  be  sub-optimal  in  the 
sense  that  an  observation/measurement-noi se  coupling  occurs,  which  can  degrade 
performance.  Homodyne-detection  interferometers  in  which  the  phase  shift  in 
one  arm  is  the  conjugate  of  that  in  the  other  arm  do  not  suffer  from  the  pre¬ 
ceding  drawback.  Overall,  however,  the  performance  of  differenced  direct- 
detection  and  homodyne-detection  interferometers  is  similar  in  single-frequency 
operation.  In  particular,  both  detection  schemes  reach  the  standard  quantum 
limit  on  position  measurement  sensitivity  in  single- frequency  interferometric 
gravity-wave  detectors  at  roughly  the  same  average  photon  number.  This  limit 
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arises  from  back  action  in  the  form  of  radiation  pressure  fluctuations  entering 
through  the  energy-phase  uncertainty  principle.  Multi -frequency  devices  can 
circumvent  this  uncertainty  principle,  as  was  illustrated  by  the  conceptual 
design  we  have  given  for  a  two-frequency  interferometer  which  can  greatly  surpass 
the  standard  quantum  limit  on  position  sensing.  This  interferometer  relies  on 
the  heterodyne  results  developed  in  [9],  and  complements  recent  work  of  Yuen  [23] 
on  the  validity  of  the  position-sensing  SQL. 
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ABSTRACT 

Degenerate  four  wave  mixing  (DFWM)  has  been  suggested  as  a  possible 
generation  scheme  for  squeezed  state  light,  A  recent  analysis  of  the 
quantum  effects  of  probe-conjugate  loss  in  backward  DFWM  has  shown  that 
such  loss  puts  an  absolute  limit  on  the  squeezing  that  can  be  obtained 
via  this  generation  scheme.  In  this  communication  we  show  that  it  is 
the  counter-propagating  beam  geometry  of  backward  DFWM  that  makes  it 
ill  suited  for  squeezed  state  generation.  On  the  other  hand,  the  nomi¬ 
nally  copropagating  beam  geometry  of  forward  DFWM  is  shown  to  alleviate 
the  absolute  probe-conjugate  loss  limit  on  squeezing. 


PACS  Headings:  42.50.  +  q,  42.65.  -  k,  05.30.  -d 


I.  Introduction 

Degenerate  four  wave  mixing  (OFWM)  has  been  suggested  by  Yuen  and  Shapiro  as 
a  possible  source  for  squeezed  state  generation^.  Their  model  was  a  simple  exten¬ 
sion  of  the  classical  description  of  DFWM  given  by  Yariv  and  Pepper^.  Quantizing 
only  the  probe  and  the  signal  beams  while  retaining  classical  descriptions  for 
the  pump  beams  and  the  nonlinear  medium,  they  showed  that  a  two-photon  coherent 
state  (TCS)  (essentially  a  minimum  uncertainty  squeezed  state)  is  obtained  by  a 
50/50  combination  of  the  phase  conjugate  reflected  beam  and  the  transmitted  probe 
beam  from  backward  DFWM.  A  recent  analysis  by  Bondurant,  Maeda,  Kumar  and 

Shapiro  has  shown  that  probe-conjugate  loss  puts  an  absolute  limit  on  the  squeez¬ 
es 

ing  that  can  be  obtained  via  backward  DFWM  .  Since  then  Reid  and  Walls  have 
given  a  fully  quantum  mechanical  treatment  of  backward  DFWM4.  Their  analysis 
neglected  the  spatial  propagation  effects  and  showed  that  pump  induced  spontane¬ 
ous  emission  limits  the  amount  of  squeezing  achievable.  In  this  communication 
we  show  that  the  absolute  limit  on  probe-conjugate  loss  is  because  the  preceding 
work  all  addressed  backward  DFWM,  which  has  a  counter-propagating  beam  geometry. 

This  geometry  is  ideal  for  correcting  phase  aberrations  via  conjugate-wave  genera¬ 
tion,  but  is  ill-suited  for  squeezed  state  generation  because  of  the  afore-mentioned 
probe-conjugate  loss  limit.  We  show  that  forward  DFWM,  which  has  a  nominally 
copropagating  nonplanar  beam  geometry,  removes  the  absolute  probe-conjugate  loss 

limit.  Such  an  interaction  geometry  has  been  applied  recently  in  studies  of  pres- 
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sure  induced  four  wave  mixing  interactions  . 

In  section  II  we  start  with  a  classical  analysis  of  forward  DFWM.  It  is 
well  known  that  large  nonlinearities  are  obtained  when  the  operating  frequency 
is  chosen  near  an  atomic  or  molecular  resonance.  Therefore,  in  section  III,  we 
develop  a  semi -classical  treatment  of  forward  DFWM  in  an  atomic  medium  con¬ 
sisting  of  an  ensemble  of  stationary  two-level  atoms.  In  section  IV  we  quantize 


the  electromagnetic  fields  and  examine  the  squeezing  behavior  of  the  output 
beams. 

II.  Classical  Equations 

Consider  the  geometry  shown  in  fig.  1.  Two  weak  waves  of  wavevectors 
k^  and  k^  propagate  at  small  angles  _+  4>/2  from  the  z  direction,  determining 
a  plane  P.  The  pump  waves  of  wavevectors  F^  and  F^  also  nominally  propagate 
along  the  z  direction;  F^  and  F^  are  obtained  from  and  Fg  by  rotating 
the  plane  P  along  the  CC'  axis.  With  this  choice  of  wavevectors  we  note  that 
F.j  +  F^  ■  F^  +  F4>  i.e.,  the  phase  matching  condition  is  satisfied. 

The  fields  are  taken  to  be  co-pol ari zed®  plane  waves  of  angular  frequency 

Ej(r,t)  *  j  Aj ( rj  )exp[i  (ut-Fj  -7)]  +•  c.c.  (1) 

~  2 

where  r.  denotes  the  distance  measured  along  k..  Following  Yariv  and  Pepper 
J  J 

we  can  derive  the  following  equations  for  coupled  modes  1  and  2 


1 ' 


(2) 


where  <  is  the  nonlinear  coupling  constant  given  by  (mks  units) 


★ 

<  =  ujx 


(3) 


A3A4/2cnQcos 


(3) 


x(3)  is  the  third  order  susceptibility  of  the  nonlinear  medium,  c  is  the 
speed  of  light  in  vacuum  and  nQ  is  the  background  refractive  index.  Equation 
(2)  has  the  following  solution 


A.(z)  =  cosh( i < | z)  A. (0)  -  i  -A-  s i nh ( ! < | z )  A,*(0),  (4a) 

i  i  i  <  i  c 

A*(z)  =  cosh(  \  <\z)  A*(0)  +  i  sinh('^'z)  A^O),  (4b) 


in  terms  of  boundary  conditions  at  z  =  0. 


-  4-  - 


III.  Coll  is  ion  less  Two  Level  Medium 

We  now  consider  an  ensemble  of  stationary  two  level  atoms  forming  the 
nonlinear  medium.  The  atoms  are  characterized  by  a  dipole  moment,  y ,  and 
longitudinal  and  transverse  relaxation  times  T-j  and  Tg  respectively.  Follow 
ing  Abrams  and  lind^  we  obtain  the  following  equations  for  the  coupled  modes 
1  and  2 


dA 


dz 

where 


zrr  =  -yai 


*  *  dA« 

•ua2  -  d r 


-y&2  +  i<  A^ 


(5) 


2a, 


uJ 


(1  +  6  ) _  J _ 

IA,V*  1**1^ 


(6) 


and 


=  2a0  J±L  U  +  D 


0  JU 

n  uj 


A  2+  A  2\2 
LA3L  +  Ia4I 


2A3A4 


As  cosU/2) 


(7) 


2  2  2 

5  =  (w-iDgjT j  is  the  normalized  detuning  from  line  center,  A$  =  ^  /T-jT^y  is 
proportional  to  the  line-center  saturation  intensity,  m  is  the  atomic  transition 
frequency,  1  =  y  AN  Tow  /2c  cfS  is  the  line-center  small-signal-field 

OO^OO 

attenuation  coefficient,  aNq  =  (N-j-N^)  is  the  equilibrium  population  difference 
in  the  absence  of  the  applied  fields,  n  is  the  saturated  refractive  index 
given  by 


and  k  is  the  magnitude  of  the  propagation  vectors  in  the  medium.  We  note  that 
the  pumps  are  nominally  copropagating,  so  that  no  spatial  averaging  along  the 
pump  direction  is  required.  Such  averaging  drastically  reduces  the  DFWM 
reflectivity  in  the  conventional  counter-propagating  pumps  geometry^. 

Equation  (5)  yields  the  following  solution. 


A-|(z)  3  e~yZ  [cosh(  |k|z)  A1  (0) 


—  sinh( | < [ z)  A*(0)] 

ii  ^ 


(9a) 


A2*(z)  3  e^2  [i  jq- sinh(  |k|z)  A^O)  +  cosh ( )  < | z)  A*(0)],  (9b) 

in  terms  of  boundary  conditions  at  z  3  0. 

IV.  Squeezed  Stace  Generation 
A.  Lossless  Case 

In  giving  a  quantum  treatment  of  backward  DFWM,  Yuen  and  Shapiro^  replace 
the  complex  field  amplitudes  A-  and  A-  with  the  photon  annihilation  and 

J  J 

creation  operators  ai  and  a.  ,  respectively,  for  j  =  1,2.  They  assume  that 
the  pump  fields  A^  and  A^  are  strong  and  hence  can  be  treated  classically. 

They  also  describe  the  medium  by  a  classical  third  order  susceptibility. 

Using  their  approach  for  our  forward  DFWM  geometry  we  replace  Eq.  (2)  with 


-  6  - 


i 


which  has  the  following  solution  for  a  z- 0  to  zs 

=L  interaction 

i 

a-j(L)  =  ya,(0)  -  iva2(0) 

(11a) 

a^L)  =  ya^O)  +  1v*a,(0). 

(lib) 

1 

Here  u  =  cosh(|<[L),  v  =  e-1 0sinh( |< JL) ,  <  s 

e10,  a,(0)  and  a2(0)  are  the 

input  field  operators  at  z  =  0.  The  outputs  at 

z  =  L  are  combined  through 

a  50/50  beam  splitter  to  generate  two  new  modes 

described  by  annihilation 

* 

operators  c  and  d  such  that 

c  =  [3,(1)  -  ia2(L)]/21/2 

(12a) 

» 

d  =  [a, (L)  +  ia2(L)]/21/2, 

(12b) 

in  terms  of  which  the  solutions  become 

w 

a 

c  =  uc.  -  vc.  +,  d  *  ud-  +  vd.  f 
in  m  »  in 

03) 

where  cin  =  [a,(0)  -  ia2(0)]/21/2 

(14a) 

a 

din  =  [8,(0)  +  ia2(0)]/21/2 

(14b) 

are  annihilation  operators  describing  field  modes  obtained  by  linear  combination 

2  2 

I  of  the  input  modes  to  the  four-wave  mixer.  Because  |yj  -|v|  =1,  it  follows 

O 

that  c  and  d  are  in  TCS  if  a,(0)  and  82(0)  are  in  coherent  states  (CS). 

B.  OFWM  in  a  lossy  medium 

I  It  was  shown  in  section  IV  A  that  modes  c  and  d  are  in  TCS.  Let  us 

j- 

c+c 

concentrate  on  the  quadrature  noise  behavior  of  mode  c.  Let  c,  *  — 2 —  and 
C2  =  -'2^—  be  the  in-phase  and  out-of-phase  quadratures  of  mode  c,  respecti vely . 

•  Then  from  Eq .  (13)  one  can  show  that  the  quadrature  variances  are 

<ac-|2>  =  <(c-)-<c1>  )2>  =  |  u-v  | 2/ 4  (15a) 

*  ^c22>  =  <(c2-<c2>)2>  =  ;.+v;2/4, 


(15b) 


when  a-j(O)  and  a^O)  are  in  CS  .  Thus  arbitrarily  large  squeezing  is  obtained 

★ 

in  c-|  when  |v|  is  made  arbitrarily  large  with  uv  real  and  positive.  Large 
values  of  |v|  have  been  shown  to  be  obtainable  in  resonant  media,  such  as 

g 

described  in  section  III  .  An  inspection  of  Eqs.  (6),  (7)  together  with  the 
defining  equations  for  y  and  v  shows  that  a  large  value  of  |v|  is  necessarily 

accompanied  by  a  large  value  of  y.  the  loss  per  unit  length  in  the  medium. 

3 

We  follow  the  approach  of  Bondurant  et  al .  to  analyze  the  effect  of  this 

probe-conjugate  loss  on  the  squeezing  obtainable  via  forward  DFWM  in  a  resonant 
medium. 

We  note  that  Eq.  (10)  can  be  obtained  from  the  effective  interaction 
Hamiltonian 

it  4. 

Hj  ■  >hv(<a1a2  +  <  a2a.|)  (16) 

using  the  Heisenberg  equations  of  motion  and  then  converting  the  temporal 
differential  equations  into  spatial  differential  equations  by  the  change  of 
variable  z  =  vt. 

In  order  to  account  for  probe-conjugate  loss  quantum  mechanically  we  adjoin 
the  system  of  Eq.  (16)  to  two  reservoirs  of  loss  oscillators10  described  by 
annihilation  operators  b^ ,  for  i  =  1  to  ®  and  s  =  1,2.  The  total  effective 
interaction  Hamiltonian  can  therefore  be  written  as 


Hj1  =  hv[<a^a2  +  <  a 2Tai  1 


+  E  [a,  z  <  o  b p  +a<-  ^ 


bp  ] , 


-  8  - 


where  represents  the  coupling  between  the  modes  of  interest,  i.e.,  a-j  and 
and  the  loss  oscillator  modes.  From  Eq.  (17),  we  obtain  two  coupled 
spatial  differential  equations  for  the  slowly- varying  operators  a^  and  a2 


da,  *•  j. 

~dT  =  "  yal  '  iK  a2  +  M2) 

dap  +  t 

“dF  =  ‘  Ya2  +  1  Kal  +  G2(z) 

where  y  is  the  loss  per  unit  length  and  Gs(z),  s  =  1 ,2  are  Langevin  noise 
operators  obeying 

yas(z)  =  v‘2  E  |<^j2  /  as(z')exp[i(u»^-u)(z'-z)/v]  dz* ,  s=l  ,2  - 
and 


(18a) 

(18b) 


(19a) 


Gs(z)  =  -  ^  E  (0)exp[i(w-aj^)z/v];  s=l  ,2 

respectively.  These  noise  operators,  under  the  Wigner-Weisskopf  approximation, 
obey  the  commutation  rule 

[Gs(z)f  Gs+,(z')]  =  2y  & ss,  6(z-z')  for  s ,  s' =1 ,2. 

The  set  of  equations  (18)  can  be  integrated  with  the  result 

a] (L)  =  e"yL[ya1(0)  -  iva2+(0)]  +  r] 

a2+(l)  =  e"YL[iva-|  (0)  +  ua2+(0)]  +  iy 

where 


(19b) 


(20) 


(21a) 

(21b) 


r  =  /  e_Y(L"z'  )[COSh(  |k  |  (L-z'))G-,  (z '  )-ie-iGsinh(  |<  ]  (L-z ' )  )G„J'(  z'  )]dz '  (22a) 

i  0  1  41 

r„  =  /  e'v^L'2  ^[-ie'lGsinh(  |<  |  (L-z'))G,  '  (z1 )  +  cosh  ( |  <  j  ( L-z' )  )G„(  z 1 )  ]dz ' .  (22b) 

0 


The  loss  per  unit  length,  >,  appearing  in  Eq .  (19a)  is  numerically  the  same 
as  that  in  Eq .  (6)  for  a  medium  consisting  of  stationary  two  level  atoms. 
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Operator  equation  (21)  reduces  to  the  classical  equation  (5)  when  expectation 
values  are  taken. 

To  calculate  the  effect  of  probe-conjugate  loss  on  squeezing,  we  construct 
new  modes  as  in  Eq.  (12)  and  evaluate  the  quadrature  variances.  As  an  example, 


c  *  e_YL[,iC.n  -  vc^n]  +  (r1  -  ir2)/21/2  (23) 

and 

<AC-j 2>  *  |2yL  (u-v)2  +  (rll+r22^2>  ’  (24) 

where  =  (r1+r1+)/2,  r22  *  (r2-r2+)/2i  and  we  have  chosen  pump  phases  such 
that  0*0.  After  substituting  Eqs.  (22a)  and  (22b)  into  Eq.  (24)  and 
evaluating  the  appropriate  moments,  we  get 


<ac-| 


and 


> 


y(2N+1] 

1 

4(y+  < 

IT 

s  y(2N+1) 

My-MT 


(25a) 


(25b) 


where  N  measures  the  initial  excitation  of  the  reservoir  modes,  i.e., 

N  *  <b^  +  (0)  b^  (0)>,  and  is  assumed  to  be  the  same  for  all  the  modes. 
Several  cases  of  interest  can  now  be  considered. 


i)  >*0,  i.e.,  the  zero  probe-conjugate  loss  limit,  in  which  Eq.  (25)  reduces  to 

Eq.  (15)  and  ideal  squeezing  is  obtained. 

1  2  2 

ii)  y^O  ,  y<  |  <  I  and  L  >>  j  - 1  p  In  this  limit  <ac2  >  ■+•  «  and  <ac^  >  ■* 

,  i.e.,  for  a  given  N  ideal  squeezing  can  be  obtained  by  making  | < |  >>  y  . 
Though  it  should  be  noted  that  pump  induced  spontaneous  emission  noise  will  limit 
this  squeezing  as  is  the  case  in  backward  DFWM  as  shown  in  Ref.  4. 


iii)  y>|<!  and  L  >>  In  this  limit  <aci2>  and  <ac?2> 


4(y+!<  i ) 


y( 2N+1 

) 

4(y-  < 

nr 

.  For  N=0,  i.e.,  when  the  loss  oscillators  are  initially  unexcited, 


we  get  <fa>  =  -fa 


1/4  and  <ac22>  =  -  j  for  y>>|< 


H) 


This  result  is  expected  here  because  loss  totally  dominates  the  nonlinear 
coupling  and  any  quadrature  noise  asymmetry  caused  by  the  latter  is  swamped 

by  the  fluctuations  introduced  by  the  former.  In  the  case  of  y=M*  a  squeezing 

7  2  2 

factor  of  2  is  still  obtained  in  <ac^  >.  The  uncertainty  product  <ac1  ><ac2  >  * 

(2N+1  )^/16(1-  '^-)>  -yg-  implying  that  a  squeezed  state  which  is  not  a  minimum 

Y  ★ 

uncertainty  state  is  generated.  Also,  since  for  our  choice  of  0,  yv  is  real, 
this  state  is  not  a  TCS  either.^ 

V.  Discussion 

The  results  of  the  previous  section  show  that  the  DFWM  beam  geometry 
plays  an  important  role  in  determining  the  squeezing  that  can  be  obtained  in 
a  realistic  experiment.  The  nonlinear  coupling  introduces  quantum  noise 
asymmetry  between  the  quadratures  of  the  interacting  modes.  In  DFVJM  this 
asymmetry  is  between  the  quadratures  of  two  different  modes  and  mode  mixing  at 
the  output  of  the  DFWM  interaction  is  required  to  obtain  new  modes,  whose  two 

quadratures  show  this  asyunetry.  Probe-conjugate  loss,  on  the  other  hand,  Intro¬ 
duces  independent  fluctuations  into  the  two  quadratures  which  are  coupled  via  the 
non-linear  interaction,  thus  tending  to  equalize  the  observed  output  beam  quadra¬ 
ture  fluctuations.  The  latter  is  the  result  of  case  (iii)  in  section  IV,  where 
loss  dominates  the  nonlinear  coupling. 

In  the  counterpropagating  geometry  of  backward  DFWM,  the  interaction 
at  any  point  couples  forward  and  backward  going  waves.  Because  of  loss, 
each  of  these  waves  has  suffered  the  noise-symmetrizing  effect  noted  above. 

It  is  the  combination  of  the  loss  with  the  non-local  nature  (forward/backward 
wave  coupling  at  all  points  in  the  interaction  medium)  that  is  responsible, 
we  believe,  for  the  severe  loss-limit  on  backward  DFWM  squeezed-state 
generation.  On  the  other  hand,  in  forward  DFWM  only  forward  going  waves 
are  coupled.  Although  loss  injects  a  symmetric  noise  contribution  at  each 
point  in  the  interaction  medium,  the  nonlinear  interaction  from  that  point 


to  the  end  of  the  interaction  squeezes  that  noise  contribution.  Thus,  with 
the  forward  interaction  and  gain  coefficient  in  excess  of  the  loss  coefficient, 
the  only  fundamental  limit  on  achievable  forward  DFWM  squeezing  will  be  due  to 
pump  induced  spontaneous  emission.  Indeed,  our  view  of  the  physics  of  this 
problem  is  supported  by  Yuen's  loss  analysis  for  DPA  squeezed  state  genera- 

O 

tion  ,  which  shows  that  in  that  forward-going  three-wave  interaction  arbitrary 
squeezing  is  obtained  for  any  y  and  N  so  long  as  |k|  can  be  made  arbitrarily 
larger  than  y. 

In  summary,  forward  DFWM  appears  more  promising  than  backward  DFWM  as  a 
squeezed  state  generator.  It  is  a  phase-matched  interaction  with  no  funda¬ 
mental  limit  on  squeezing  due  to  probe-conjugate  loss.  Moreover,  we  expect 
there  will  be  differences  in  the  limits  on  obtainable  squeezing  set  by  pump- 
induced  spontaneous  emission  in  forward  and  backward  DFWM,  because  of  the 
different  physics  of  their  spatial  propagation  characteristics,  as  described 
above.  Furthermore,  since  all  the  beams  are  propagating  in  roughly  the  same 
direction,  the  interaction  is  not  velocity  selective.  The  participation  of 
all  velocity  groups  results  in  a  very  large  nonlinear  interaction.  Both 

experimental  and  theoretical  investigations  of  forward  DFWM  have  recently 
12 

been  published.  Phase  conjugation  and  sub-Doppler  resolution  due  to  strong 
saturation  have  been  reported. 
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O.F,  Walls,  This  work  was  supported  by  the  Office  of  Naval  Research  Contract 
N0001 4-81 - K-0662 . 


References 

1.  H.P.  Yuen  and  J.H.  Shapiro,  Opt.  Lett.  £,  334  (1979). 

2.  A.  Yariv  and  D.M.  Pepper,  Opt.  Lett.  J_,  16  (1977). 

3.  R.S.  Bondurant,  M.  Maeda,  P.  Kumar  and  J.H.  Shapiro  in  "Coherence 
and  Quantum  Optics,  V,"  eds.  L.  Mandel  and  E.  Wolf,  pp.  767-774, 
Plenum  Press,  NY,  1983;  R.S.  Bondurant,  P.  Kumar,  J.H.  Shapiro 
and  M.  Maeda,  submitted  to  Phys.  Rev.  A. 

4.  M.D.  Reid  and  D.F.  Walls,  submitted  to  Optics  Communications. 

5.  Y.  Prior,  Appl .  Opt.  J9_,  1741  (1980);  J.A.  Shirley,  R.J.  Hall,  and 
A.C.  Eckbreth,  Opt.  Lett.  5 ,  380  (1980);  A,R.  Bogdan,  Y.  Prior, 
and  N.  Bloembergen,  Opt.  Lett.  6^,  82  0981 L;  Y.  Prior,  A.R.  Bogdan, 
M.  Dagenais,  and  N.  Bloembergen,  Phys.  Rev.  Lett.  4£,  111  (19.81). 

6.  We  assume  the  polarization  to  be  perpendicular  to  plane  P  and  ignore 
the  complicating  effects  due  to  a  small  angle  $  between  the  electric 
field  vectors  of  the  pumps  and  iji/2  between  the  electric  field  vector 
of  the  pump  and  the  probe. 

7.  R.L.  Abrams  and  R.C.  Lind,  Opt.  Lett.  2,  94,  (1978);  3,  205  (1978). 

8.  H.P.  Yuen,  Phys.  Rev.  An,  2226  (1976). 

9.  R.C.  Lind  and  D.G.  Steel,  Opt.  Lett.  6.,  554  (1981). 

10.  W.H.  Louisell,  "Quantum  Statistical  Properties  of  Radiation,"  Wiley, 
New  York,  1973. 

11.  H.P.  Yuen,  Phys.  Lett.  56A,  105  (1976). 

12.  Y.  Prior  and  E.  Yarkoni,  Phys.  Rev.  A  28,  3689  (1983);  A.  Khyzniak, 
V.  Kondilenko,  Yu.  Kucherov,  S.  Lensik,  S.  Odoulov,  and  M.  Soskin, 

J.  Opt.  Soc.  Am.  A  1_,  169  (1984). 


Forward  DFWM  geometry. 


APPENDIX  II 


QUANTUM  NOISE  AND  EXCESS  NOISE  IN  OPTICAL 
HOMODYNE  AND  HETERODYNE  RECEIVERS 


Jeffrey  H.  Shapiro 


Department  of  Electrical  Engineering  and  Computer  Science 
and  Research  Laboratory  of  Electronics 
Massachusetts  Institute  of  Technology 
Cambridge,  MA  02139 


Abstract 

A  parallel  development  of  the  semiclassical  and  quantum  statistics  of 
multi -spatiotemporal  mode  direct,  homodyne,  and  heterodyne  detection  using  an 
ideal  (except  for  its  sub-unity  quantum  efficiency)  photon  detector  is  presented. 
Particular  emphasis  is  placed  on  the  latter  two  coherent  detection  configurations. 
The  primary  intent  is  to  delineate  the  semiclassical  theory's  regime  of  validity 
and  to  show,  within  this  regime  of  validity,  how  the  quantum  theory's  signal 
quantum  noise,  local  oscillator  quantum  noise,  the  quantum  noise  incurred 
because  of  sub-unity  detector  quantum  efficiency,  plus  (for  heterodyning  only) 
image  band  quantum  noise  produce  the  quantitative  equivalent  of  the  semiclassical 
theory's  local  oscillator  shot  noise.  The  effects  of  classical  fluctuations 
on  the  local  oscillator,  and  the  recently  suggested  dual-detector  arrangement 
for  suppressing  these  fluctuations,  are  treated.  It  is  shown  that  previous 
studies  of  this  arrangement  have  neglected  a  potentially  significant  noise 
contribution . 
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I .  INTRODUCTION 

In  coherent  optical  detection  [1]  -  [3],  the  optical  field  to  be  measured 
is  combined  on  the  surface  of  a  photodetector  with  the  field  of  a  strong 
local  oscillator  laser  whose  center  frequency  is  offset  by  an  amount  Av  from 
that  of  the  signal  field.  The  detection  scheme  is  referred  to  as  optical 
homodyning  if  Av=0,  and  optical  heterodyning  if  Av  *vjp>0,  with  v^p  being  the 
intermediate  frequency  in  the  latter  case.  For  both  schemes,  electrical 
filtering  of  the  photocurrent  is  used  to  select  the  beat  frequency  components 
in  the  vicinity  of  Av,  yielding  an  output  that  contains  a  frequency  translated 
replica  of  the  signal  field  components  that  were  coherent  in  space  and  time 
with  the  local  oscillator  field.  Heterodyne  detection  is  now  widely  employed 
in  coherent  COg  laser  radars  [4],  [5],  and  is  being  vigorously  researched  for  use 
with  semiconductor  injection  lasers  in  fiber  optics  [6]  -  [8]  and  space 
communications  [9],  [10].  Performance  analyses  in  these  areas  routinely  employ 
the  semiclassical  statistical  model  for  photodetection,  which  implies  that  the 
fluctuations  observed  in  coherent  optical  detection  with  signal  and  local 
oscillator  fields  of  perfect  amplitude  and  frequency  stability  comprise  an  additive 
white  Gaussian  noise,  representing  local  oscillator  shot  noise. 

It  has  long  been  known  [11]  that  the  semiclassical  statistics  for  photodetec¬ 
tion  are  quantum  mechanically  correct  only  when  the  total  field  illuminating  the 
detector  is  in  a  Glauber  coherent  state  or  a  classically  random  mixture  of  such 
states.  Inasmuch  as  ordinary  light  sources,  including  lasers  and  light  emitting 
diodes,  obey  this  classical  state  condition,  there  is  no  need  to  abandon  the 
semiclassical  approach  in  the  vast  majority  of  photodetection  sensitivity 
calculations.  However,  non-classical  light  has  been  generated  via  resonance 
fluorescence,  as  confirmed  by  observations  of  its  photon  anti-bunchi ng  [12]  and 
sub-Poissonian  behavior  [13]  in  direct  detection.  Moreover,  there  is  great 
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theoretical  interest  in  squeezed  states  {also  called  two-photon  coherent  states) 

[14],  [15],  which  are  non-classical  states  of  considerable  potential  for  optical 
communications  [16]  -  [19]  and  precision  measurements  [20]  -  [23],  For  these 
states,  the  quantum  theory  of  photodetection  is  essential,  and  coherent  optical 
detection  schemes  are  the  most  interesting. 

In  [18],  Yuen  and  Shapiro  developed  the  quantum  descriptions  of  single¬ 
detector  optical  homodyne  and  heterodyne  receivers.  They  employed  a  quasimono- 
chromatic  approximation,  and  assumed  a  coherent  state  local  oscillator, 
corresponding  to  perfect  local  oscillator  amplitude  and  frequency  stability. 

Within  these  limitations,  complete  statistics  for  mul ti-spatiotemporal  mode  detection 
are  available  from  [18].  More  recently  [22],  the  fact  that  photodetectors  respond 
to  photon  flux  rather  than  power  [24]  has  been  used  to  relax  somewhat  the 
quasimonochromatic  approximation  in  [18].  Thus,  were  high  power,  highly  stable 
local  oscillator  lasers  available  at  all  wavelengths  of  interest,  the  quantum 
photodetection  theory  of  [18]  would  provide  a  sufficiently  general  foundation  for 
all  optical  homodyne  and  heterodyne  sensitivity  calculations.  Unfortunately , 
such  is  not  the  case. 

Driven  by  heterodyne-detection  problems  arising  from  the  excess  noise  of 
semiconductor  injection  lasers,  Yuen  and  Chan  [25]  proposed  a  dual -detector 
arrangement  for  coherent  optical  detection,  akin  to  the  balanced  mixer  concept  of 
microwave  technology  [26],  [27].  They  gave  a  direct  quantum  analysis  of  single- 
mode  dual-detector  homodyning,  showing  that  local  oscillator  quantum  and  excess 
noises  can  be  balanced  out,  hence  alleviating  injection  laser  problems  that  would 
have  plagued  a  single-detector  system.  In  subsequent  work  by  Chan  and  his 
collaborators ,  the  basic  dual -detector  excess  noise  cancellation  concept  was 
demonstrated  experimentally  [28],  and  a  variety  of  non-ideal  device  effects  (quantum 
efficiency  mismatch,  etc.)  were  analyzed  using  semiclassical  mul ti -temporal  mode 
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techniques  [29].  Also,  Schumaker  [30]  has  shown  that  the  dual-detector  single- 
mode  homodyne  arrangement  is  better  than  single-detector  homodyning  for  making 
non-classical  squeezed  state  observations,  as  a  result  of  its  ability  to  cancel 
out  local  oscillator  excess  noise. 

Because  the  quantum  treatments  in  [25]  and  [30]  are  confined  to  single-mode 
situations,  and  the  multi-mode  results  in  [28],  [29]  are  in  essence  semiclassical , 
there  is  as  yet  no  fully  quantum  treatment  of  multi -mode  dual -detector  coherent 
optical  reception.  This  paper  will  develop  such  a  model  by  generalizing  the 
results  of  [18].  Simple  explicit  representations  for  all  of  the  relevant  output 
terms  in  coherent  optical  detection  with  a  strong  but  classically  random  local 
oscillator  field  will  be  derived.  It  will  be  seen  that  the  previous  dual¬ 
detector  analyses  [25],  [28]  -  [30]  neglect  excess-noise  modulation  of  the  signal 
and  quantum  noise  terms,  and  the  first  of  these  modulation  effects  may  significantly 
degrade  output  signal-to-noise  ratio  in  some  circumstances.  Moreover,  because 
of  the  calculational  power  afforded  by  [18]?  our  rather  general  quantum  results  are 
more  directly  comparable  with  those  of  the  multi-mode  semiclassical  theory  than 
are  the  more  limited  results  of  [25],  [30].  Indeed,  that  comparison  is  the 
primary  purpose  of  this  paper. 

The  paper's  core,  Section  II,  is  a  parallel  development  of  the  semiclassical 
and  quantum  statistics  of  mul ti-spatiotemporal  mode  direct,  homodyne,  and  heterodyne 
detection  using  an  ideal  (except  for  its  sub-unity  quantum  efficiency)  photon 
detector.  The  formulation  therein  for  the  coherent  optical  detection  schemes  will 
assume  perfectly  stable  local  oscillators  in  the  semiclassical  models,  and  the 
corresponding  coherent  state  local  oscillators  in  the  quantum  models.  We  use 
Section  II  to  delineate  the  semiclassical  theory's  regime  of  validity,  and  to  show, 
within  this  regime,  how  the  combination  of  the  quantum  theory's  signal  quantum 
noise,  local  oscillator  quantum  noise,  the  quantum  noise  incurred  because  of 
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sub-unity  detector  quantum  efficiency,  plus  (for  heterodyning  only)  image  band 

quantum  noise  produce  the  quantitative  equivalent  of  the  semiclassical 

theory's  local  oscillator  shot  noise.  In  Section  III  we  address  coherent 

optical  detection  with  classically  random  local  oscillators.  The  technique  of 

iterated  expectation  is  used  to  readily  obtain  both  semiclassical  and  quantum 

results  for  this  case.  Single-detector  and  dual-detector  systems  are  considered, 

and  our  results  are  compared,  in  the  case  of  dual-detector  quantum  homodyning  to 

those  of  [25],  [30].  Finally,  in  Section  IV  we  briefly  discuss  the  implications 

of  our  work  for  squeezed  state  generation  experiments,  which  is  the  application 

that  motivated  our  analysis. 

II.  SEMICLASSICAL  VS.  QUANTUM  PH0T00ETECTI0N 

The  central  element  of  all  the  photodetection  configurations  we  will  consider 

is  shown  in  Fig.  1.  It  is  a  surface  photoemitter  with  active  region 

x  *  (x,y)  z  in  the  z=0  plane,  illuminated  by  a  quasimonochronomatic  (center 

frequency  v  )  paraxial  scalar  electromagnetic  wave  from  the  half  space  z<0  over 

an  observation  time  interval  t  e  T.  This  detector  is  assumed  to  have  a  constant 

quantum  efficiency  n  ver  the  frequency  band  containing  the  illuminating 

field.  The  output  of  the  detector  is  a  scalar  current  density  J(7,t)  for 

x-  £  A  ,  t  s  T.  As  will  be  described  below,  the  field  characterization  we 
d 

must  employ  for  the  illumination  is  either  classical  or  quantum  mechanical,  according 
to  whether  semiclassical  or  quantum  photodetection  statistics  are  sought.  Although 
we  shall  neglect  internal  time  constant  and  noise  effects,  which  are  present  in 
real  detectors,  our  direct  detection  results  will  be  applicable  to  photomultiplier 
tubes  (for  which  the  current  gain  permits  internal  noise  to  be  overcome)  at 
post-detection  bandwidths  up  to  the  reciprocal  anode  response  time  of  the  tube. 
Furthermore,  our  results  will  be  applicable  to  coherent  optical  detection  systems 
using  semiconductor  photodiodes  (for  which  the  mixing  gain  overcomes  the  internal 
noise)  up  to  the  post-detection  bandwidth  of  the  detector.  No  particular  loss  of 
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generality  is  entailed  by  the  use  of  scalar  rather  than  vector  fields,  with 
the  caveat  that  all  the  coherent  optical  detection  work  herein  presumes  that 
the  actual  signal  and  local  oscillator  fields  are  co-polarized.  Finally,  by 
appropriate  spatial  integrations,  we  can  collapse  our  current  density  observation 
to  photocurrent  observations  for  a  single  detector  or  a  multiple-detector 
array. 

A.  Direct  Detection 

In  direct  detection,  the  electromagnetic  field  to  be  measured  comprises  the 
entire  illumination,  and  the  basic  observation  quantity  is  the  current  density  J(7,t). 
Semi  cl  as  si  cal  Model  Let  E^(x,t)  be  the  positive-frequency  complex 

field  (V/m  units)  associated  with  the  classical  scalar  electric  field  incident 
on  the  detector,  i.e.,  E^(7,t)  is  the  analytic  signal  of  this  electric  field. 
8ecause  of  our  quasi  monochromatic  assumption,  the  Fourier  transform  1  of  E^ 

E(+)(7,v)  =  [dt  E(+)(7,t)  ej2irvt  (1) 

> 

is  non-zero  only  for  |v  -  vQ|  <_  B,  where  the  bandwidth  B  is  much  less  than  the 
center  frequency  Because  of  our  paraxial  assumption,  the  short  time  average 

power  density  falling  on  the  point  7  at  time  t  is 

1(7, t)  =  (ce0/2)  E(')(7,t)E(+)(7,t),  (2) 

( _ \  f  +1  * 

where  c  is  the  speed  of  light,  £Q  the  permittivity  of  free  space,  and  E  i(E  ') 

★ 

is  the  negative- frequency  complex  field,  with  denoting  complex  conjugate. 

"he  standard  semiclassica 1  photodetection  model  [31],  in  our  notation, 
presumes  that  J(7,t)  is  a  conditional  space-time  Poisson  impulse  train  with 
rate  function  u(x,t)  =  RI(7,t)/e  where  e  is  the  electron  charge,  and  ? 

is  the  detector's  responsivity  (A/W  units)  at  the  illumination's  center 

* 

frequency  .■  .  This  means  that: 
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1)  the  current  density,  which  is  of  the  form 

J(x,t)  *  J[e5(x-xn)6(t-tn)t  (3) 

n 

has  shot  effect  noise,  i.e.,  it  consists  of  instantaneous 
emissions  of  an  electron  charge  e  at  the  random  space-time 
points  {(7n,tn)  :  x^A^,  tneT}; 

2)  conditioned  on  knowledge  of  the  rate  function  (u(x,t)  :  xeA^.teT}, 
the  number  of.  photoemissions  occurring  within  a  spatial  region 

A'C  Ad  during  a  time  interval  V  C_  T  is  a  Poisson  random  variable 

r  _ .  r  _ 

with  mean  value  j  dx  ’  |  dt  u(x,t); 

A'  T’ 

3)  conditioned  on  knowledge  of  {u(7,t)  :  xeA^  ,  teT} ,  the  photoemissions 
occurring  in  disjoint  spaial  regions  A',  A"CAd  are  statistically 
independent  processes, 

4)  conditioned  on  knowledge  of  {u(x\t):  jTeA^.ter},  the  ohotoemissions 
occurring  in  disjoint  time  intervals  T'  ,T"0’  are  statistically 
independent  processes. 

Even  though  the  semiclassical  theory  of  photodetection  employs  classical 

fields,  it  is  customary  to  recognize  in  this  theory  that  light  of  frequency  v 

is  quantized  into  photons  of  energy  hv  ,  where  h  is  Planck's  constant.  Thus, 

for  the  quasimonochromatic  case  at  hand,  the  responsivity  is  ordinarily  written 

as  R  =  en/hv  ,  in  terms  of  the  detector's  quantum  efficiency  n  and  the 
o 

photon  energy  at  the  field's  center  frequency,  so  that  u(x\t)  =  nI(7,t)/hvQ. 
In  fact,  because  we  are  concerned  with  detectors  that,  quantum  mechanically, 
respond  to  photon-flux  density  rather  than  power  density  [24],  [22],  it  is  more 
proper  to  write 

u(x,t)  -  nl ( x , t ) , 


(4) 
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where  I  h  is  the  classical  photon-flux  density 


Iph(x.t)  =  E  (x,t)E(x,t) 


obtained  from  the  (photons/s)  /<L/m  units  positive- frequency  complex  field 


E(x,t)  =  [  dv(ce  /2hv)  1/2  E(+)(7,v)e“j2,rvt. 


For  all  practical  purposes  in  the  semi  classical  theory,  with  quasi  monochromatic 
light  we  can  use  u(7,t)  =  nI(7,t)/hvQ  and  y(7,t)  =  nlpfl(7,t)  inter¬ 
changeably.  This  amounts  to  using  v  =  vQ  in  the  square-root  term  of  (6), 
an  approximation  whose  validity  is  guaranteed  by  (1).  In  the  quantum  theory, 
even  with  quasi  monochromatic  light,  it  is  critical  to  employ  the  photon-flux 
density  formulation,  see  [22]. 

Quantum  Model  In  the  quantum  photodetection  theory,  the  classical  positive- 

frequency  complex  field  E^(7,t)  is  replaced  by  a  positive-frequency  field 

operator  -  E^(x,t),  whose  quantum  state  is  specified  by  a  density  operator  o. 

The  quasi  monochromatic  and  paraxial  conditions  of  the  semiclassical  theory 

become  conditions  on  the  density  operator,  namely,  that  the  excited  (non- vacuum 

state)  modes  of  E^(x,t)  lie  at  frequencies  within  B  of  v  and  propagate 

at  small  angles  to  the  z  axis.  As  in  [17],  [18],  [22],  we  shall  regard  the 

current  density  J(7,t)  as  a  classical  quantity,  corresponding  to  the  macroscopic 
2 

output  of  the  quantum  measurement  performed  by  the  detector  of  Fig.  1  on 
the  field  E^  '(x,t).  To  provide  an  explicit  representation  of  this  quantum 
measurement,  we  must  first  develop  the  quantum  effective  photon-flux  density. 

Let  us  convert  E^(x,t)  to  a  photon-units  field  operator  by  defining 
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(cf.  Eqs.  (1),  (6)  ) 


E(+)(x,v)  *  dt  E^(x,t)  , 

E(x,t)  *  |  dv  (ccQ/2hv) E(+)(7,v)e~j2lTvt  . 


Equation  (8)  defines  the  same  basic  field  operator  used  in  [16]  -  [18],  [22]. 
We  can  made  a  modal  expansion 


E(x,t)  =  Ian5n(x,t),  xeA.,  teT  (9) 

n 

of  this  operator,  where  {an>  are  modal  annihilation  operators  satisfying  the 
commutation  rules 


*  »•  CvO  ■ 


and  (5n>  are  a  complete  orthonormal  set  of  classical  functions  over 
xeA^,  teT.  In  Eq.  (10),  the  fap+}  are  the  adjoints  of  the  {an>;  they 
are  modal  creation  operators. 


For  a  detector  of  sub-unity  quantum  efficiency  we  must  adjoin  to  (9)  a  fictitious 


field 


2vac^x,t^  =  -  Cn^n^x’^^’  Xe^d’  ^e^"’ 


where  (c)  are  modal  annihilation  operators  that  commute  with  {a  } 
n  n 

and  {an+},  viz. 
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tvs3  =  [viJ  *  =0- 


cs-O 


nm. 


(12) 


The  fields  E  and  Evac  are  quantum-mechani  cally  independent,  with  the  latter 
having  all  its  modes  in  the  vacuum  state.  In  terms  of  E  and  Ew._ 

the  effective  photon-flux  density  operator  for  the  detector  is 


Iph(x,t)  =  E,+(x,t)E' (x,t)  , 


(13) 


with 

E* (X.t)  =  n1/2E(7,t)  +  (1-n)1/2Evac(7,t) 


(14) 


The  representation  theorem  of  quantum  photodetection  [18,  theorem  1]  can 
now  be  stated  (in  our  notation)  as  follows.  The  classical  current  density 
J (X.t)  obtained  from  photoemissive  detection  measures  the  quantum  operator 

J(x,t)  =  eiph(x,t)  ,  (15) 

i.e.,  it  is  proportional  to  the  effective  photon-flux  density.  In  somewhat 
more  detail  this  means  (cf.  the  semiclassical  case): 

1)  the  current  density  obeys 

J(x,t)  .  I  e  S(x-J )  ,  (16) 

n 

so  it  is  still  a  collection  of  instantaneous  emissions  of  an  electron 
charge  at  random  space-time  points  ((x  ,t  )}  ; 
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2)  if  F(J(7,t))  is  an  arbitrary  functional  of  the  current  density,  then 
the  classical  average  of  this  random  variable  <F(J(7,t))>  equals 
the  quantum  average  tr(p ' F(J(7,t))) ,  where  tr  denotes  trace  and 

o'  ’p®pu>,  gives  the  joint  density  operator  for  E  and  E„_,  in 
terms  of  the  density  operator  p  for  E  and  the  vacuum-state 

pM,  *  ®  1 0x0 1  density  operator  for  E 
n 

Note  that  we  cannot  dispense  with  the  vacuum  state  field  E..__  unless  n*!, 

VflC 

A 

even  though  its  average  value  obeys  tr(pw„  Ew  (x,t))  *  0  regardless  of  the 
value  of  n  •  This  is  because  the  zero-point  fluctuations  (vacuum-state  quantum 
noise)  in  Eyac  can  contribute  to  F(J(7,t)).  Indeed  the  noise  in  J(7,t)  has 
nothing  to  do  with  the  snot  effect  associated  with  the  discreteness  of  the 
electron  charge.  Rather,  it  is  the  quantum  noise  in  E'  being  observed  through 
measurement  of  the  effective  photon-flux  density. 

A 

Comparison  Let  us  suppose  that  the  density  operator  for  E  is  a  classical  state 
i  .e. , 


(17) 


wi  th 

a_ |o>  3  a  |a>  (18) 

n  ,  n  - 

defining  the  multi -mode  Glauber  coherent  states  of  the  field  E  in  terms  of 

♦ 

the  modal  expansion  (9),  and  P(a;o  )  being  a  classical  probability  density 
function  P(a;a  )_>0,  jd  a  P(a;a  )  =1.  It  was  shown  in  [18]  that  this  is  a 

necessary  and  sufficient  condition  for  the  semiclassical  statistics  to  be 
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quantitati vely  correct.  In  particular,  under  this  condition  the  quantum  model 
predicts  that  J(x,t)  is  a  conditional  space-time  Poisson  impulse  train  with 
conditional  rate  function 

u(x,t)  »  nE*(x,t)E(x,t)  (19) 


for 

E(x,t)  =  l  anCn(x,t)  *  <ajE(x,t)|a>  ,  (20) 

n  ' 

the  average  field  illuminating  the  detector  given  the  state  of  E  is  the  multi-mode 
coherent  state  |a  >  . 

To  illustrate  the  above  behavior,  let  us  examine  the  statistics  of  the 
a 

observed  photon  count 


N 


dx  dt  J(x,t) 
Ad  Jr 


(21) 


assuming  single-mode  illumination,  and  n  *  1.  In  the  semiclassical  theory  we 
shall  take 


E(x,t)  3  [a/(AdT)1/2]e'j2,rvot,  ^Ad,t£T  , 


(22) 


where  a  is  a  complex-valued  random  variable  with  probability  density  function 
p(a),  Ad  is  the  area  of  Ad,  and  T  is  the  duration  of  T.  We  then  obtain 
Handel's  rule  [32]  for  the  probability  distribution  of  N 

=  [d2a  p (a)  ( |a|2n/n! )  exp  ( - | a | 2)  , 


Pr[N  3  n] 


(23) 
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viz.  N  is  a  conditionally  Poisson  random  variable.  Equation  (23)  gives 
the  mean  and  variance  of  the  observed  ohotocount  to  be 


<N> 


d 2a  p  (a)  |a  | 2 


(24) 


and 


var(N)  =  <N>  +  var  (|a|2),  (25) 

respectively,  where  the  first  term  on  the  right  in  (25)  represents  shot 
noise  and  the  second  term  on  the  right  in  (25)  represents  excess  noise. 

In  the  quantum  theory  we  let  C-j ( x*. t )  ■  (AdT)"1^2e"j2irvot  be  the  mly 
excited  mode  in  (9),  so  that  the  density  operator  for  E  is 

°3oln®|0><°l  (26) 

for  o-|  the  density  operator  of  mode  1.  We  then  find  for  the  probability 
distribution  of  N  [11],  [17] 

Pr[N  *  n]  =  <n|p1  |n>  ,  (27) 

where 

*  0.* 

a.| ' a.|  i n>  =  n  | n>  (28) 

defines  the  photon  number  states  of  the  first  mode  of  E.  If  o ^  is  the 
classical  state. 

£>1  =  jd2a  p(a)  | a><a |  ,  (29) 

with  p(a)  being  the  probability  density  from  the  semiclassical  theory 


_T_  ^  •  JL.  -* 


(cf.  Eq.  (17)),  then  (27)  reduces  to  (23)  as  expected.  Thus,  in  this  case  the 
semi  classical  theory  is  quantitatively  correct  in  its  prediction  of  the  photon 
counting  probability  distribution.  It  is  nevertheless  physically  incorrect,  in 
that  it  ascribes  the  photon  counting  fluctuations  to  shot  noise,  whereas  they  are 
actually  a  manifestation  of  the  illumination  field's  quantum  noise.  For 
example,  were  *  { k><k |  where  |k>  is  the  k-photon  number  state  (a 

non-classical  state),  then  we  would  get 

Pr[N=n]  =  5nk  (30) 

from  (27),  whence 

<N>  *  k  ,  (31) 

and 

var(N)  =  0  ,  (32) 

for  the  photon  count  mean  and  variance.  Here  the  field  state  is  an  eigenket  of 
our  observation  operator,  so  there  is  no  uncertainty  in  the  measurement 
outcome.  This  sub-Poissonian  behavior  cannot  be  obtained  from  the  semiclassical 
theory,  because  for  all  p(a)  the  excess  noise  term  in  (25)  will  be  non-negative, 
forcing  var(N)>.  <N>  to  prevail. 

3.  Homodyne  (Detection 

The  configuration  v*.  shall  consider  for  single-detector  multi -spatiotemporal 
mode  homodyne  detection  is  shown  in  Fig.  2.  The  signal  field  to  be  detected  is 
combined,  through  a  lossless  beam  splitter  of  intensity  transmission  £,  with  a 
perfectly  stable  local  oscillator  field  on  the  surface  of  the  Fig.  1  photodetector. 
The  resulting  current  density,  J^om(x\t),  is  our  homodyne  detection  output, 
whose  statistics  we  shall  characterize  below.  By  spatial  integration  of  our  results 
over  the  detector's  active  region  A  , 


we  can  use  our  model  to  describe  single- 


detector  homodyning;  the  extension  to  dual-detector  homodyning  will  be  made  in 
Section  III. 


Semi  classical  Model  The  total  classical  photon-units  complex  field  incident  on  the 

e 

photodetector  is  given  by 

E(xft)  *  £1/2E$(7,t)  +  (l-e)1/2ELQ(7,t)  .  (33) 

for  7eAd  ,  teT,  in  terms  of  a  (potentially  random)  weak  signal  field 

Es(7,t),  and  a  deterministic  strong  local  oscillator  field  ELg(7,t).  The 
latter  has  a  classical  photon-flux  density 

*phl0^x’^  3  ^LO  ^x,^^E^g(x,t)  ,  (34) 

that  greatly  exceeds  that  of  the  former 

Iphs^»t^  *  Es#(7.t)Es(7,t)  ,  (35) 

for  xeAd,teT  .  Thus,  the  rate  function  driving  the  photodetector  is,  from 
(4).  (5),  (33)-(35),  approximately 

u(7,t)  *  n[(l-c)IphL0(7,t)+2[e(l-£)]1/2Re(Es(7,t)EL0*(7,t))]  .  (36) 

It  then  follows,  from  the  Central  Limit  Theorem  for  high  density  shot  noise  [33], 
[34],  that  at  very  large  values  of  the  local  oscillator  classical  photon  number 


NphL0 


dx 
A 


dt  I 
r 


phLO 


(x,t) 


(37) 


the  homodyne  detection  current  density  ^orr/*’^  is  a  conditional  Gaussian 
process.  Specifically,  conditioned  on  knowledge  of  the  signal  field 
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{ Eg ( x , t )  :  xeA^,  tcT}  ,  J^om(x«t)  is  the  sum  of  three  current  densities: 

1)  a  homodyne-mixing  current  density  signal  term  2en[e(l-e)]1/2Re(Es(x,t)ELO*(x\t) ) 

2)  a  direct-detection  local  oscillator  bias  current  density  en(l-e)IphLgCx,t)  ; 


3)  a  local  oscillator  shot  noise  current  density,  which  is  a  zero-mean 
spatiotemporal  non-stationary  white  Gaussian  noise  process  ^Shot^x,1:^ 
covariance  function 


<Jshot(xl  ,t:l)Jshot(x2’V: 


'n(  1  -e ) IphL0^xi  ^ (xi ‘X2)<S( - 1 


In  order  to  connect  the  preceding  mul ti-spatiotemporal  mode  formulation  with 
more  familiar  single-detector  multi-temporal  mode  results,  let  us  consider  the 
statistics  of  the  single-detector  homodyne  photocurrent 


“  !  dx  JhnJX> t)  , 


assuming  that 


EL0(x’t}  =  (PL0/hvoAd)V2  e'j2,Vot  -  (40) 

corresponding  to  a  normally-incident  plane  wave  local  oscillator  of  power 
Here  we  find  that,  conditioned  on  knowledge  of  the  signal  field,  i^Qm(t) 
comprises  a  signal  current 


1  s i q ( t )  =  2en[PL0£p-s)/hvoAd]1/2Re(  j  dx  E$(  x,t)eJ'  2™o*)  , 


plus  a  zero- frequency  bias  current 


bias 


en(  1  ■ 


•€)pL0/hvo 


(42) 


plus  a  zero-mean  stationary  white  Gaussian  process  shot-noise  current 
i$hot(t)  with  sPectral  density  eib-as  (A2/Hz)  . 

The  signal  current  is  a  frequency-transl ated  (to  baseband)  replica  of  the 
normally- incident  plane  wave  component  of  E^Cx.t)  that  is  in  phase  with  the 
local  oscillator  field.  The  bias  current  is  the  zero-frequency  photocurrent 
produced  by  the  local  oscillator  field.  The  noise  current  is  the  local -oscillator 
shot  noise,  whose  spectrum  follows  the  well  known  Schottky  formula  [35]. 
quantum  Model  In  the  quantum  model,  £q.  (33)  becomes  an  operator-valued  expression 

E(x,t)  *  £1/2Es(7,t)  +  (l-e)1/2ELQ(7,t),  (43) 

giving  the  field  operator  E  that  drives  the  detector  in  terms  of  the  signal 
field  operator  E<.  and  the  local  oscillator  field  operator  Elq  .  The  density 
operator  o  for  E  is  assumed  to  be 

3S  ®  °L0  (4a) 

where  o,  is  an  arbitrary  signal  field  density  operator  and  -,n  3  a  ><a  [ 

5  LU  '10  -L0 

is  a  multi-mode  coherent  state  local  oscillator  density  operator.  The  latter 

corresponds  to  a  mean  local  oscillator  field 


Eiq(x,t)  5  <-L0^L0^*’t^--l0>  =  -  1L0n  • 


when  E^g  is  expanded  using  the  mode  set  •  -:n }  as  was  done  for  E 


(45) 

in  Eq .  (9) . 


The  strong  local  oscillator  condition  of  the  quantum  theory, 


tr(ps£sf(x,t)Es(x,t))  «  |ELQ(x,t)  |2 


is  assumed  to  prevail  (cf.  Eqs.  (34),  (35)),  with  a  very  large  average  local 
oscillator  photon  number 


dx  [  dt  | E.  n(x,t)  |2»  1  , 


Ad  T 

(cf.  Eq.  (37)). 

To  obtain  the  effective  photon-flux  density  operator  measured  by  the  detector 
we  adjoin  to  E  from  (43)  a  quantum-mechanically  independent  vacuum-state 
field  operator  i  „  .see  Eqs.  (11)  -  (14).  We  can  now  give  a  fully  quantum 
characterization  of  the  classical  homodyne  current  density  J.  (7,t)  ,by 
translating  the  results  of  [18,  theorem  Z]  into  our  notation.  The  strong  local 
oscillator  condition  implies  that  this  classical  current  density  measures  the 
quantum  operator 

Jhom(x,t)  s  en(l-e)EL0(x,t)EL0(x,t)  +  2e[n(T-e )]1/2Re  {[(ne)1/2Es(7,t) 

+  (l-n)1/2Evac(7,t)]  •  EL(J+(7,t)  }  .  (48) 


Moreover,  because  N^g>>! ,  the  local  oscillator  direct  detection  term  in  (48), 

en(l-s  )ELg(7,t)ELg(7,t)  ,  yields  classical  observation  values  comprising  a  bias 

current  density  en(l-e) |E^g(x,t) |  plus  a  local  oscillator  quantum  noise  current 
density,  which  is  a  zero-mean  spatiotemporal  non-stationary  white  Gaussian  noise 
process  ^LOq^’^  w1  covariance  function 

<JlWxl’VJL0<^x2,V> 


[er(l-c)]  EL0(xrt,)  -  • '  t.--; 
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Furthermore,  under  this  same  condition^,  the  second  term  on  the  right  in  (48) 

1  /  9  ~  _  ★  _ 

simplifies  to  a  homodyne-mixing  signal  operator  2en[e(l-e)]  '  Re(E<.(x,t)E^g(x,t)) 

plus  a  sub-unity  quantum  efficiency  (n<l)  quantum-noise  current  density.  The  latter 
current  density  is  a  zero-mean  spatiotemporal  non-stationary  white  Gaussian 
noise  process  J  (7,t)  with  covariance  function 

VuC 

<Jvac(xl  = 

e2n(l-n)(l-e)|EL0(71,t1)|25(71-72)6(t1-t2)  ;  (50) 

is  statistically  independent  of  .  ^1US*  the  classical  homodyne  current 

density  J.  m(7,t)  measures  the  operator 

Jhom^x,t^  =  en(l-e) (ELg(x,t) |  +  JL0q^x,t^ 

+  2en[£(l-e)]1/2  Re(Es(7,t)ELg(7,t) )  +  Jya(;(7,t).  (51) 

The  first  term  on  the  right  in  (5T)  is  the  local  oscillator  bias  current  density, 

the  second  term  is  the  classical  representation  of  the  local  oscillator's  quantum 

noise  contributed  by  the  en  (i-e)Eio  E^g  measurement,  and  the  last  term  is  the 

classical  representation  of  the  n<l  quantum  noise  contributed  by  the 

2e[-( l-n)(l-c ) ]^2Re(Eyac(7, t)ELg(7,t) )  measurement.  The  signal  field 

contribution  to  J.  (7,t)  cannot  be  simplified  further  without  knowledge  of 
nom 

the  density  operator  .  In  general,  this  term  will  contribute  signal  field 
quantum  noise  to  the  homodyne  observation,  as  will  be  seen  below. 

Comparison  To  facilitate  comparison  of  the  semiclassical  and  the  quantum  theories 
of  homodyning,  we  shall  restrict  our  consideration  to  the  single  detector  case. 

First,  we  need  the  quantum  characterization  of  the  homodyne  photocurrent  (39), 
which  can  be  obtained  by  spatial  integration  of  the  results  just  presented.  We 
assume  a  normally-inci dent  plane  wave  mean  local  oscillator  field 


and  we  find  the  measures  the  operator 


^hom^  ~  Si  as  +  Soq^^ 

i  2  t 

+•  2en[PLge(1-e)/hvoAd]'l,/^Re(  dx£s(x,t)  e  0  )  +  i 


vac 


(t) 


(53) 


Here,  ts.  is  given  by  (42),  and  i,_  and  i  „  are  statistically 

OiaS  LUh  vaC 

independent  zero-mean  stationary  white  Gaussian  noise  processes  with  spectral 
densities  en(l-e  ) i 5-,- as  and  e^"n^bias  •resPect'i vely  • 

Physically,  1’bl-as  is  the  local  oscillator  bias  current,  iLQq  is  the  local 
osci llator  quantum-noise  current,  and  ivac  is  the  n<l  quantum-noise  current. 
Equation  (53)  differs  from  the  semiclassical  description  in  two  respects:  the 
homodyne-mixing  signal  term  involves  the  quantum  field  operator  E<.  rather 
than  the  classical  field  E^  ;  the  noise  in  the  homodyne  observation  is  a 
combination  of  local  oscillator  quantum  noise,  n<l  quantum  noise,  and  signal 
quantum  noise,  rather  than  simply  being  local  oscillator  shot  noise.  We  know, 
from  the  di  rect  detection  discussion,  that  the  semiclassical  photodetection  model 
is  quanti tati vely  correct  if  the  density  operator  0  for  the  field  E  illuminating 
the  detector  represents  a  classical  state.  This  situation  occurs,  under  (44), 
if  and  only  if  j«-  .the  signal  field  density  operator,  is  a  classical  state 


-s 


'  2  * 

=  |d  .S'^S  )  I  =ts><rotS  ■ 


(54) 


for  the  multi-mode  signal  field  coherent  state  in  modal  expansion  of 

E^  similar  to  Eq.  (9),  with  P^  being  a  classical  probability  density.  When  (54) 

applies,  the  homodyne-mixing  signal  term  in  (53)  can  be  given  a  classical 

representation  akin  to  that  employed  for  the  E  mixing  term  in  going  from 

vac 


(48)  to  (51).  In  particular,  for  a  classical  r.ignal  field  state,  the  quantum 
theory  of  homodyning  predicts  that 


where 

Es(7,t)  =  <as|Es(7,t) |as>  (56) 

is  the  classical  mean  signal  field  when  the  state  of  is  |<x->  .  and 
ij  (t)  is  a  zero-mean  stationary  white  Gaussian  noise  current  of  spectral 
density  enei^s  that  is  statistically  independent  of  i^g  and  i'vac 

The  classical  field  is,  in  general,  a  random  process  with  probability 

★ 

density  Ps(as;as  )  in  modal  expansion  form.  The  current  iSq(t)  is  the 
classical  representation  of  the  coherent  state  signal  field  quantum  noise  as 
observed  through  the  measurement  operator  (53).  Note  that 


i,(t)  '  W'1  *  V'1  +  V'1 


is  a  zero-mean  stationary  white  Gaussian  noise  process  of  spectral  density 
e i 5 i as  ’in  quantitative  agreement  with  the  semiclassical  is^Qt(t)  result. 

Of  course,  the  interpretation  of  the  origin  of  the  noise  in  homodyning  is 
different  in  these  two  theories.  Local  oscillator  shot  noise  is  a  semiclassical 
fiction;  the  noise  seen  in  homodyne  detection  (with  an  ideal  local  oscillator) 
is  local  oscillator  quantum  noise,  plus  n<l  quantum  noise,  plus  signal  quantum 


-22- 


noise.  Moreover,  in  the  limit  ne  -*■  1  with  n(l-e)NLg  >>  1  .the  two  former 
contributions  disappear,  and  homodyning  gives  a  direct  quantum  measurement  of 
the  signal  field  component  that  is  coherent  in  space  and  in  phase  with  the 
local  oscillator  [18].  It  is  this  characteristic  that  makes  homodyning 
attractive  for  squeezed  state  applications  [16]  -  [19]. 

C.  Heterodyne  Detection 

The  configuration  for  single-detector  mul ti-spatiotemporal  mode  heterodyne 
detection,  shown  in  Fig.  3,  mimics  that  employed  for  homodyne  detection.  The 
only  differences  are  that  the  signal  field  is  centered  at  frequency  vQ  +  vjp  , 
the  local  oscillator  is  centered  at  frequency  \>Q  ,and  passband  filtering  of 
the  current  density  is  used  to  select  beat  frequency  components  in  the  vicinity 
of  the  IF  frequency  Vjp  (\jjpT>>l  will  be  assumed).  The  bandwidth  B  of 
the  signal  field  will  be  taken  to  be  much  less  than  v^p  .and  we  shall  concern 
ourselves  with  characterizing  the  statistics  of  the  current  density 
Jhet(x\t).  The  results  we  need  are  easily  developed  by  injecting  the  frequency 
offset  vjp  into  the  preceding  homodyne  work. 

Semiclassical  Model  In  Eq.  (33)  let  us  make  the  frequency  offset  of  the  signal 
field  explicit  by  writing 

E$(x,t)  =  Fs(7,t)e'j2:T(vo+'JIF)t  ,  (58) 

where  F^  is  a  baseband  complex  signal  field  of  bandwidth  B  .  The  results 
following  (37)  now  provide  the  semiclassical  statistics  for  heterodyning,  namely, 
conditioned  on  knowledge  of  the  baseband  signal  field  (F^(7,t)  : 

7eA^ ,  teT)  ,  Jy.  ^(7,t)  is  the  sum  of  three  current  densities: 

1)  a  heterodyne-mixing  current  density  signal  term 

2e-1[£(l-O]1/2Fe(Fs(7,t)e‘j2:r(vo+'JIF)tEL0*(7,t))  ; 
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2)  a  direct-detection  local  oscillator  bias  current  density 
en(l-s)IphL0(x,t)  ;  and 

3)  a  local  oscillator  shot  noise  current  density  characterized 


by  (38). 

The  single-detector  heterodyne  photocurrent 

d*  Jhet(x.t)  , 


W(t>  = 


(59) 


a 

assuming  is  given  by  (40),  then  comprises  a  signal  current 


isig(t)=2en[PL0e(l-£)/  hvQAd]1/2Re(|  dxFs(x,t)e'j2ir',IFt) 


(60) 


plus  a  bias  current  1.  .  from  (42),  plus  a  zero-mean  stationary  white  Gaussian 
process  shot-noise  current  1shQt(t)  with  spectral  density  eih. 


bias* 


The  heterodyne  current  (59)  is  thus  a  frequency  translated  (from 
Jo  ^IF  to  VIF^  version  of  the  normally  incident  plane  wave  component  of 
plus  the  usual  bias  and  shot  noise  terms,  Because  of  the  frequency  offset 
'Jrp  between  the  signal  and  the  local  oscillator  fields,  both  the  in-phase  and 
quadrature  (relative  to  the  local  oscillator)  components  of  the  signal  field 
contribute  to  the  output  observations. 

Quantum  Model  Here  we  suppose  that  the  only  non-vacuum  state  modes  of  the 

field  operator  E^  lie  within  a  bandwidth  B  of  the  frequency 

However,  because  of  zero-point  fluctuations,  the  quantum  version  of  (58)  is 


E,(7,t)  =  Mx,t)e-j2T(VvIF>t  +  FT(I,t)e-j27r(W)t:  . 


(61) 


where  F^  and  F^  are  baseband  complex  signal  and  image  field  operators. 
Physically,  the  image  band,  being  v^p  Hz  below  the  local  oscillator’s  frequency, 
contributes  quantum  noise  to  even  when  it  is  unexcited  [18],  [22], 

We  shall  assume  that  F^  and  Fj  are  quantum-mechanical ly  independent,  with 
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f-i 

* 


the  latter  having  all  its  modes  in  the  vacuum  stated  Wc  now  find,  from  . 
the  quantum  homodyning  work,  that  J^e£(7,t)  measures  the  operator 

Jhet(x,t)  =  en(l-€)|EL0(7,t)|2  +  JL(Jq(7,t) 

+  2en[e(l-£)]1/2Re(Fs(7,t)e"j2Tr(uo^IF)tEL0*(7,t)) 

+2en[E(l-e)]1/2Re(FI(7,t)e‘j2ir(VvIF)tEL0*(7,t)J  +  Jvac(7,t)  ,  (62) 


where  E.  n,  J.  n_  ,  and  J  are  as  given  in  (51).  We  can  use  the  vacuum-state 
LU  LU^  vac 

nature  of  F^  to  obtain  the  classical  representation 


2en[£(l-£)]1/2Re(FI(7,t)e‘j2ir(VvIF)tEL0*(7,t)) 

-  JIq(7,t)/21/2  ,  (63) 

where  Jjq  is  a  zero-mean  non-stationary  white  Gaussian  classical  process 
corresponding  to  the  image-band  quantum  noise,  with  covariance  function 


<JIq(xl »ti)Jiq(x2*t2^>  s 

(en)  c  ( 1  -e  )  |  E^q( , t^ )  |  3 ( x-j -*2 ) ®  ( t-j -tj ) 


(64) 


Thus,  the  quantum  description  of  the  si ngle-detector  heterodyne  photocurrent 
i^  (t)  from  (59)  is  that  it  measures  the  operator 


WU>  ■  'bias  *  Wc>  +  iIq‘t)/2'/Z  *  \ac(t) 

-  2e-[PLOc(l-€)/hvQAd]1/2Re(  [  d7F$(7,t)e‘j2'TVIFt) 


(65) 


where  ibias  •  1  LOq  •  and  1yac(t)  are  aS  in  and  1  is  a  zer0 

mean  white  Gaussian  noise  process  (the  classical  representation  of  image-band 

quantum  noise)  of  spectral  height  ene 7bi as ’ 

Comparison  Suppose  the  density  operator  for  F^  is  a  classical  state,  i.e., 

its  density  operator  pF  obeys 

S 

p  =  |d2as  PF  ^aS;aS  ^  las><aS^  * 

‘5  j  -  S''  " 

where  PF  is  a  classical  probability  density,  and  |cu>  is  the  multi-mode 
S 

Glauber  coherent  state  for  the  modal  expansion 


F-(x,t)  *  i a  5  (x.t)«JMV',IF)t  ,  (67) 

b  n  Sn 

with  l  denoting  summation  over  modes  Cn  lieing  within  bandwidth  B  of 
n 

frequency  v0+vrp.  Here  we  can  obtain  a  classical  representation  of  the 
F  term  in  (65)  which  reduces  the  quantum  description  of  the  heterodyne 
photocurrent  to 

=  'bias  +  'LOq^^  +  +  W*1 


+2enLPLQe(l-€)/hv0Ad]1/2Re(  Jdx  F$(x,t)e  _jZlTVIFt) 


+  iSq(t)/2,/4  , 

.  ,,1/2  ,  .  /?l/2  .  , 

where  the  total  noise  current,  iLOq  +  iIq/*  +  ^q7  ' 
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is  a  zero-mean  white  Gaussian  process  with  spectral  density  ei.  .  , 

0 1  dS 

in  quantitative  agreement  with  the  semi  classical  theory,  and 

Fs(7,t)  =  <as|Fs(7,t)  |as>  (69) 

is  the  classical  baseband  signal  field  envelope  F^  associates  with  the 
coherent  state  |o^>  .  Note  that  half  of  the  quantum  noise  entering 

’het  comes  through  the  signal  field  operator  F^  and  the  other  half 

-  Q 

comes  through  the  image  field  operator  F. . 

III.  EXCESS  NOISE  EFFECTS  AND  DUAL-DETECTOR  OPERATION 

In  this  section  we  shall  extend  the  results  of  Section  II  for  coherent  optical 
reception  to  include  classical  excess  noise  on  the  local  oscillator  field  and 
dual-detector  operation.  It  is  convenient  to  begin  with  a  presentation  of 
dual -detector  results  in  the  absence  of  excess  noise. 

A.  Dual -Detector  Coherent  Optical  Reception 

Suppose  the  homodyne/ heterodyne  configurations  of  Figs.  2  and  3  are  augmented 
by  the  use  of  another  quantum  efficiency  n  detector  on  the  previously  unused 
output  port  of  their  beam  splitters,  see  Fig.  4.  We  take  the  classical  output 
field  for  this  port  to  be 

E(7,t)  =  -(l-e)1/2Es(7,t)  +  £1/2ELQ(7,t)  (70) 

in  the  semiclassical  model,  and  use  the  corresponding  operator-valued 
expression  in  the  quantum  model.  Rather  than  treat  the  full  mul ti -spatiotemporal 
mode  situation,  we  shall  restrict  our  attention  to  the  photocurrents 
i.|(t)  and  i^(t)  obtained  by  spatial  integration  of  the  current  densities 
J-|(x,t)  and  J2(7,t)  produced  by  detectors  1  and  2.  We  shall  assume  a 
perfectly  stable  (i.e,  deterministic)  classical  local  oscillator  field 
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EL0(x>t)  “  FL0(x,t)e-j2irV  (71) 

with  baseband  complex  envelope  in  the  semiclassical  model,  and  a 

Glauber  coherent  state  quantum  local  oscillator  with  mean  field 

tr(pL0EL0(x,t))  =  FLQ(x,t)e"^27rvot  (72) 

with  baseband  complex  envelope  FLg  in  the  quantum  model.  Under  these 
conditions  the  results  of  Section  II  can  be  used  to  show  that  the  following 
statistics  apply. 

Homodyne  Detection  In  homodyning,  the  signal  field  is  centered  on  v  ,  so, 
because  of  (71), (72),  it  is  convenient  to  introduce  baseband  signal  complex 
envelopes  via 

Es(x,t)  3  Fs(7,t)e“j2,cV  (73) 

and 

Es(x,t)  =  Fs(7,t)e_j2lrVot  ,  (74) 

for  the  semiclassical  and  quantum  cases,  respectively.  Now  we  have,  semi- 
classically,  that 

i-i(t)  3  en(l-fi)  j  dx|F.  n(7,t)  |2 


+  Zen[e(l-e)]"^2Re(  j  dxF«-(x,t)  F.  -  (x,t)) 
■Ad 
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and 


for  the  homodyne  photocurrents,  where  nsh0tl^  and  nshot2^  are  s'ta't'is't1cally 
independent  identically  distributed  zero-mean  stationary  white  Gaussian  noise 
processes  of  unity  spectral  density.  Equations  (75)  and  (76)  have  the  usual 
bias  plus  mixing  signal  plus  local-oscillator  shot  noise  interpretation.  Note 
that  the  beam  splitter  phase  shift  between  the  output  ports  forces  the  mixing 
signals  to  be  180°  out  of  phase.  Also,  the  independence  of  the  local - 
oscillator  shot  noises  follows  because  they  are  generated  from. deterministic 
illumination  of  two  different  detectors. 

For  the  quantum  case,  we  have  that  i^(t)  and  i^(t)  measure  the  operators 

i-j(t)  »en(l-«)j  dx|FLQ(7,t)  |2 
Ad 

+  2en[e(l-e)]1/2Re(j  dxFs(x,t)FLQ*(7,t)) 

Ad 

+  en(l-e)[j  d7[FL0(7,t)i2]1/2  nLOq(t) 

Ad  f 

+  e[n(l-n)(l-e)  dx j F^ x, t) | 2]] /2nvacl ( t)  ,  (77) 

Ad 

i7(t)  3  enc !  dx[F  n(x,t) |2 

;Ad 

-2enCe(l-e)]1/2Re(J  d7  Fs(x,t)FLQ*(7ft)  ) 


and 


+e[n(l-n)£  dx|F  0(x,t) I  ]  '  nvac2(t) ,  (78) 

Ad 

where  nvacl^’  nvac2^  are  stat’st‘'ca!^>'  independent  identically 

distributed  zero-mean  stationary  unity- spectrum  Gaussian  processes.  The  familiar 
bias  plus  quantum  mixing  signal  plus  loci! -oscillator  quantum  noise  plus 
n<l  quantum  noise  interpretation  applies  to  (77),  (78).  As  in  the  semiclassical 
model,  the  mixing  term  appears  180°  out  of  phase  in  the  two  photocurrents.  No 
such  phase  shift  appears  on  the  n^  term,  as  this  noise  arises  out  of  the 
direct  detection  of  ELq.  Indeed,  except  for  scale  factors,  the  local-oscillator 
quantum  noise  contributions  to  i-j(t)  and  ^(t)  are  completely  correlated.  The 
n<l  quantum  noises  are,  on  the  other  hand,  statistically  independent  because  they 
arise  from  different  detectors.  Finally,  when  the  signal  field  is  in  a  classical 
state  these  quantum  results  can  be  shown  to  be  in  quantitative  agreement  with 
the  foregoing  semiclassical  formulas. 

Heterodyne  Detection  For  heterodyning  we  use  (58),  rather  than  (73),  to  introduce 
a  baseband  signal  complex  envelope  for  the  semiclassical  analysis.  We  then  find 
that 

i -j (t)  ■  en(l-e)  j  dx|FLQ(x,t)  |2 
‘Ad 

+2en[£(l-c)]1/2Re(  |  d7Fs(7,t)FL0*(7,t)e'j2’rvIFt) 

jAd 

+«Cn(l-e)f  d7|F.n(7,t)|2]1/V.ot1(t)  , 


(79) 


(30) 


i2(t)  *  ene  dx|FQ(x,t)| 

a 

-2en[e(l-e)J1/2Re(j  d7F$(7,t)FL0*(7,t)e"j2*vIFt) 


+e[ne  dx|FL0(x,t)|  ]  ^nShot2^^  ’ 
Ad 


with  interpretations  as  given  following  (75),  (76).  In  the  quantum  case  we 

use  (61)  instead  of  (74)  and  obtain  the  measurement  operators 

i-j  ( t)  »  en(l-e)  I  dx|  FL0(x,t)  | 2 
Ad 

+2en[e(1-e)]1/2Re(j^  d7Fs(7.t)FLJ(7,t)e'j2™IFt) 
d 

+en(l-e)  [f  dx|FLQ(x.t) |2j1/2  nLQ  (t)  . 

Ad 

+e[n(l-n)(l-e)[  dx| FLQ(7,t) |2]1/2nvacl (t) 

d 

+en[£(l-€)2'1jA  d7|FLQ(I,t)|2]1/2nIq(t)  ,  (81) 


for  detector  1  and 

i 2 ( t )  =  fine |  dxjFL0(x,t) |2 

d 

-2en[e(l-£)]1/2Re(|  dxF.(x,t)F  *(7,t)e"j2irvIFt) 
’Ad 

■^fine [  j  dxjFL0(x,t)  |2l1/2nL0q(t) 

'Ad 

+fi[n( 1 -n)e  j  d7|  FLQ(7,t) j 2]1/2nvac2(t) 
-«n[e(l-r)2_1f  d7!F,n(7,t)l2]1/2nT  (t)  , 


(82) 


for  detector  2.  In  Eqs  (81),  (82)  the  Interpretations  and  coinnents  following 
(77),  (78)  are  applicable.  The  noise  n^(t),  which  represents  image-band 
quantum  noise,  is  another  zero-mean  stationary  unity-spectrum  white  Gaussian 
process.  It  is  statistically  independent  of  n^(t)  and  nvacj^  f°r 
j  =  1,2,  and  appears  with  a  sign  reversal  in  i^  and  ^  because  it  arises 
from  the  mixing  term  involving  Fp 

B.  Local -Oscillator  Excess  Noise 

The  extension  of  the  results  of  Section  II I A  to  incorporate  classical 
excess  noise  on  the  local  oscillator  is  extraordinarily  simple,  because  of  the 
form  the  preceding  results  have  been  cast  in.  Specifically,  for  the  semi  classical 
theory  we  need  only  make  the  baseband  local  oscillator  complex  envelope  F^g 
in  (71)  a  complex- valued  random  process  with  known  statistics.  Then  the  homodyne 
and  heterodyne  results  of  the  semiclassical  theory,  namely  Eos.  (75),  (76) 
and  Eqs.  (79),  (80),  respectively,  become  conditional  statistics  assuming  FL0 

Q 

is  known.  Unconditional  statistics  follow,  via  iterated  expectation  [38], 
from  averaging  over  the  local  oscillator  fluctuations,  as  will  be  illustrated 
below.  In  a  similar  manner,  classical  local -oscillator  excess  noise  can  be 
injected  into  the  quantum  model  by  making  o^g  a  classical-state  density  operator 
for  which  F^g,  the  average  baseband  local -oscillator  complex  envelope  given  the 
local  oscillator  is  known  to  be  in  the  multi-mode  coherent  state  |a^g>,  1S  3 
complex- valued  classical  random  process.  The  quantum  homodyne  and  heterodyne 
results,  Eqs.  (77),  (78),  and  (81),  (82),  respecti vely,  are  now  conditional 
characterizations  given  F^g.  Unconditional  statistics  are  again  obtained  by 
averaging  over  the  local  oscillator  fluctuations.^ 

To  illustrate  our  excess  noise  results,  and  compare  them  with  relevant  prior 
work  [25],  [28]- C 30 ] ,  we  shall  consider  a  single  spatial  mode/multi -temporal 
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mode  local  oscillator,  for  which  F^q  in  the  semi  classical  theory  and  F^g  in 
the  quantum  theory  are  both  of  the  form  (pLg(t)/hvoAd)  '  exp(-j$^g(t))  ,  where 

PLQ(t)  and  iLg(t)  are  classical  random  power  and  phase  fluctuations.  For 
convenience,  we  shall  assume  that  these  fluctuations  are  the  polar  decomposition 
of  a  stationary  comp! ex-Gauss i an  random  process.  We  shall  also  assume  that  the 
signal  field,  in  both  the  semiclassical  and  quantum  pictures,  is  statistically 
independent  of  the  local  oscillator.  Finally,  we  shall  limit  our  consideration 
to  the  differenced  output  currents  i ^ ( t )  -  i2(t). 

Homodyne  Detection  Under  the  preceding  conditions  we  have  the  semiclassical 
result 


i 1 ( t )  -  i 2( t)  3  en(l-2£)PL0(t)/hvo 

+4en[£(l-€)PL0(t)/hvoAd]1/2Re(|  d7Fs(7,t)eJ\o(t)) 


+e[nH-£)PL0(t)/hvo]  /2nshotl(t) 

-e[nePLg(t)/hvQ]  nShot2^  ’ 


and  the  quantum  result 


i -j  C t )  -  i £ ( t )  =  en(l-2e)PL0(t)/hvo 

+4er,rs(l-t)PL0(t)/hvo  Ad]1/2Re(j  d7F$(7,t)eJ\o^ ) 


+en(l-2£)(PL0(t)/hvo)l/SL0q(t) 


+e[n( 1 -n) ( 1 -e )P^g( t)/hvQ]  nvaci 


-eCn(l-n)ePL0(t)/hvo]  nvac2(t) 
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In  both  (83)  and  (84),  the  first  term  on  the  right  equals  a  mean  bias  current 
en(l-2e)<PLQ(t)>/hvQ  plus  a  local  oscillator  power- fluctuation  excess 
noise  en(l-2e)(PLQ(t)-<PLQ(t)>)/hvQ.  8oth  of  these  are  exactly  nulled  when 

the  beam  splitter  is  50/50,  i.e.,  when  t~  1/2.  The  second  term  on  the  right  in 
(83)  and  (84)  is  the  homodyne-mixing  signal  current;  local  oscillator  randomness 
both  amplitude  and  phase  modulates  this  term.  The  remaining  terms  in  the 
semiclassical  result  (83)  are  the  shot  noises,  now  modulated  by  local  oscillator 
power  fluctuations.  The  remaining  terms  in  the  quantum  result  are  the  local 
oscillator  quantum  noise  and  the  n<l  quantum  noises;  these  too  are  modulated 
by  the  local  oscillator  power  fluctuations.  Note  that  when  s  =  1/2  the  local 
oscillator  quantum  noise  contribution  vanishes. 

Let  us  further  specialize  the  quantum  results  by  supposing  that  the  only 

•  -1/2 
excited  mode  of  F<.  is  the  monochromatic  plane-wave  pulse  (A^T)  for 

7EAd  ,  teT  ,  and  that  a$  is  the  annihilation  operator  for  this  mode.  Matched 

filtering  of  the  differenced  output  currents  then  yields  a  measurement  of 

•T 

M  s  e'1 jQ  [i i ( t)  -  i 2 ( t ) Jdt  ,  (85) 

where  normalization  by  the  electron  charge  has  been  used,  for  convenience,  to 
make  the  observation  values  dimensionless.  We  assume  that  the  mean  function 
and  covariance  function  of  the  stationary  complex-Gaussian  local-oscillator 
random  process 

l( t)  3  (PL0U)/hvoAd)1/2  exp(-j*LQ(t))  (86) 

are 

»  [0--')<NLO>/AdT]1/2  ,  (37) 


and 
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K^(t)  =  (y<NL0>/AdT)k(T) 


(88) 


respectively,  where  is  the  average  number  of  local  oscillator  photons 

1/2 

present  over  AdxT  ,  y  is  the  fractional  root-mean-square  (rms)  local 
oscillator  amplitude  fluctuation  level,  and  k(r)  is  a  real -valued 
normalized  covariance  (k(0)  =1).  It  then  follows  that 

<M>  =  n(l-2£)<N^Q> 

+  4n[e(l-e)(l-y)<NL0>J1/r2  <a$1>  ,  (89) 

and 


<AM2>  =•  n(l-n)<N.  n>  +  n2(l-2e)2<Min>  +  16n2e(l-e)(l-y)<N,  n><Aa  2> 


10 


10 


10  uaSl 


+  [n(l-2e)<NL0>/T] 


fT 


.2, .2 


k  (r)  +  2y(l-y)k(t)]  (T-|t| 


+  ( 2n/T) 2 ( N2e )  y[e  ( 1  -e ) <N^Q> 3 ( 1  - y)  ] 1 7 2  2  <a$1  >  jVk  ( t)  (T- 1  r  | ) 


+  (2n/T)2s(l-e)y<N^g>[T2+2<a2  a^> 


-T 

dTk(r)(T-|r|)]  , 


(90) 


-T 


give  the  mean  and  variance  of  the  M  measurement  from  which  a  signal-to-noise 
ratio 


SNR; 


M 


<M>  /<aM  > 


(91) 


may  be  calculated.  In  Eq.  (89),  the  first  term  is  the  average  local  oscillator 
bias  current  contribution,  and  the  second  term  is  the  average  signal  field 
mixing  term  contribution,  where  a$1  =  (as+as*)/2  for  a$  the  annihilation 
operator  of  the  sole  excited  Fs  mode.  In  Eq.  (90),  the  first  term  is  due  to 
the  n<l  noises  nvac]  and  nvac2  >  the  second  term  is  due  to  the  local 
oscillator  quantum  noise  nLQ  ,  the  third  term  is  the  signal  field  quantum 
noise,  the  fourth  term  is  due  to  the  local -oscillator  power  fluctuations,  and 


the  last  terms  are  due  to  the  random  modulation  of  the  mixing  current  by  the 
local  oscillator  fluctuations. 

The  previous  dual -detector  homodyne  studies  of  Yuen  and  Chan  [25]  and 
Schumaker  [30]  assume  e  *  1/2,  ( 1 -y)<N^q>>>1 ,  and  a  slowly  fluctuating  local 

oscillator  (corresponding,  in  our  case,  to  k(t)  *  1  for|T|<T).  In  this  limit  both 
prior  studies  find  (in  our  notation) 


SNR* 


<*51* 


<aas-|  >+(l-n)/4n 


(92) 


whereas  we  obtain 


SNR*m 


-  ? 


3 

<  4  a  ^  >+(  1  -n  )/4n  ( 1  -y  )+y(  1  -^a^a^ )  /4(  1  -y ) 


(93) 


For  small  fractional  rms  local -oscillator  amplitude  fluctuations  (y«1), 

(93'  differs  from  (921  because  of  an  additional  noise  term  in  the  denominator 
that  is  approximately  y(1  +■  2<a$  as>)/4.  Physically,  this  term  arises  from  the 
random  local -oscillator  modulation  of  the  mixing  current,  an  effect  neglected  by 
the  earlier  studies.  In  order  for  this  term  to  be  insignificant  compared  to  the 
signal  quantum  noise  of  a  c  .ierent  state  (<aasl2>  =  1/4),  we  require  that 

4, 

Y<as  as>  <<1  .  (94) 

i.e.,  the  fractional  rms  local  oscillator  amplitude  fluctuation  must  be  much 

smaller  than  the  square  root  of  the  reciprocal  of  the  average  number  of  signal 

field  photonsP  This  requirement  becomes  even  more  stringent  if  a  squeezed  state 

*  7 

is  being  probed,  for  which  <Aasl  >  <1/4  prevails. 
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In  addition  to  exhibiting  the  potentially  significant  random  modulation 


of  the  mixing  term,  our  formulation,  (84),  shows  another  effect  suppressed 
in  [25]  and  [30].  This  is  the  random  amplitude  modulation  of  the  local  oscillator 
and  n  <  1  quantum  noises  by  the  classical  amplitude  noise  of  the  local 
oscillator.  Although  this  modulation  does  not  explicitly  enter  the  signal-to- 
noise  ratio  ,  it  does  make  the  last  three  terms  in  (84)  non-Gauss i an  random 
processes,  an  effect  which  will  modify  digital  communication  error  probability 
calculations  somewhat. 

Heterodyne  Detection  The  semi  classical  description  for  the  differenced  output 
currents  in  heterodyne  detection  is 

1 T ( t )  -  i 2 ( t )  =  6n(l-2£)PLQ(t)/hv0 
+4en[£(l-£)PL0(t)/hvoAd]1/2Re(|  d7Fs(x,t)e*j(2irvIFWL0(t) ) ) 


*e[n(l-£)PL0(t)/hvoJ  nshot! ( t) 


-e[nePLQ(t)/hvo]  n^t)  , 


and  the  quantum  description  is 


i-j(t)  -  i  2  ( ^ )  =  2n(l -2e  )P^g(t)/hvQ 


MenU(l-OPLO(t)/hv0  Ad]1/2Re(j  dxFs(7,t)e'j  (2™IFt_\o(  } ) 

Ad 


^(l-2£)(PL0(t)/hvo)'^nL0q(t) 


^[-(l-n)(l-£)PL0(t)/hvo]l/2nyacl(t) 


+en[2£{1-£)PL0{t)/hvo]1/2nIq(t)  .  (96) 

These  results  differ  from  the  corresponding  homodyne  results,  (83)  and  (84),  in 
only  two  respects.  First,  the  mixing  terms  (second  terms  on  the  right  in  (95) 
and  (96))  beat  the  signal  field  to  an  intermediate  frequency  not  baseband, 
and  so  they  sense  both  quadratures  of  the  signal  field.  Second,  the  quantum 
result  (96)  gains  a  noise  contribution  from  the  image  band  quantum  noise  through 
nj  (t)  .  The  local  oscillator  excess  noise  (and  its  cancellation  when  e=  1/2) 
and  the  random  modulation  of  the  signal  and  noise  terms  by  the  local  oscillator 
fluctuations  thus  continue  to  be  present  in  the  heterodyne  case,  i.e.,  the 
interpretations  given  for  the  homodyne  situation  aoply  here  as  well.  Once  again, 
the  relevant  previous  work  on  dual  detector  systems  [25],  [28],  [29]  does  not 
include  all  the  effects  contained  in  our  treatment;  the  random  local  oscillator 
modulation  of  the  signal  and  noise  terms  is  absent  in  the  above  analyses. 

As  an  illustration  of  these  omissions,  let  us  compare  our  semiclassical  answer  (95) 
assuming  a  deterministic  monochromatic  plane-wave  pulse  signal  F<-Cx,t)  * 

-1/2  — 

3S^d"n  for  xeAq,t cT,  with  the  corresponding  e-1/2,  equal  quantum 

efficiency  result  of  Abbas  and  Chan  [29].  The  latter  claim,  in  our  notation,  that 
the  differenced  output  currents  consist  of  a  mean  current 

<(i](t)  -  i  2  ( t )  )>*  2enT'1<NL0>1/2Re(ase‘j2:r'JIFt)  ,  (97) 

embedded  in  an  additive  zero-mean  white  Gaussian  noise  process  with  bilateral 
spectral  density 

S(f)  =  e2n<NL0>/T  .  (98) 

We  have,  using  (86)  -  (87)  in  (95),  that  the  differenced  output  currents  consist 
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of  a  mean  current 

<(T1(t)  -  i2(t))>=2enT‘1L(l-Y)<NL0>]1/2Re(ase':i2irvIFt)  ,  (99) 

plus  a  conditionally  non-stationary  zero-mean  white  Gaussian  shot  noise  process 
that,  given  the  local  oscillator  power  waveform,  has  covariance  function 

K(t,s)  =  e2n(PL0(t)/hvQ)5(t-s)  ,  (100) 

plus  a  signal  dependent  zero-mean  stationary  Gaussian  noise  process 

i"(t)  =  2en(Ad/T)1/2Re[as(^(t)  -  m  )Yj2*',IFt]  ,  (10T) 

with  covariance  function 


K^„^i,(t)  =  2(en/T)2Y<N^Q>  |a^  | 2k(x)  COS(2irvjpr) 


(102) 


When  Eg. (100)  is  averaged  over  the  P^q  statistics  it  reduces  to  a  stationary 
white  noise  spectrum  (98)  ,  however  the  random  PLQ  fluctuations  make  the 
noise  non-Gaussian ,  albeit  in  a  minor  way  if  y<<l.  The  noise  current 
i ” ( t )  comes  from  the  random  modulation  of  the  mixing  term  and  may  present  a 
significant  degradation.  Consider  a  high  guality  (y<<1),  slowly  fluctuating 
local  oscillator  (k(-)  -  1  for  |r|<_T)  and  the  matched  filter  processor  generating 


■K1 


M  =  e*':  ( i 1 ( t )  -  i2(t) )21/2cos[2mvIFt- arg(as)] 


(103) 


then  the  Abbas  and  Chan  model  gives  a  signal -to-noi se  ratio 


SNRm 


(104) 


whereas  we  have  that 


-39- 


I 

E 


1 
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2n  I  <*s  1 2 

SNR^  «  - s—— -  .  (105) 

1  +•  nYla^l4, 

As  in  the  quantum  homodyne  example  given  earlier,  at  high  average  detected 
signal  levels  there  is  a  very  stringent  requirement  on  local  oscillator 
amplitude  fluctuations  if  SNR  degradation  is  to  be  avoided. 

IV.  DISCUSSION 

At  this  point,  we  have  clearly  established  how  the  quantum  theory  for 
coherent  optical  detection  subsumes  the  familiar  semiclassical  statistics  in 
a  natural  way.  We  have  also  seen  that  the  quantum  approach  is  essential  for 
studying  the  photodetection  statistics  of  non-classical  field  states.  There 
is  now  considerable  interest  in  a  particular  class  of  non-classical  states, 
called  the  two-photon  coherent  states  [14]  or  the  squeezed  states  [15].  These 
states  are  in  essence  minimum  uncertainty  product  states  for  the  quadrature 
components  of  the  photon-units  field  operator  E(x,t).  In  particular,  for  a 
single  field  mode  with  annihilation  operator  a,  the  two-photon  coherent  state 
j3;u,\>>  obeys  the  eigenket  relation 

(ua  +  vaT)|s;u,v>  *  b|3;u,v>  ,  (106) 

2  2 

where  s.u.v  are  complex  numbers  and  u.v  satisfy  |u|  -  |v|  =  1.  With 

a-j  =  (a  +  aT)/2  and  s  (a  *  a  )/2j  denoting  the  quadrature  components  of  a, 

we  then  find  that  the  state  |8;u,v>  gives 

<aa]2>  =  |  u  -  v  1 2/4 ,  (107a) 

and 

<  >  =  |u  +  v|2/4. 


7* 


(107b) 


When  u,v  are  real  valued,  (107)  implies  that  je;u,v>  satisfies  the 
Heisenberg  relation 

<  Aa]2>  <Aa22>  >  1/16  (108) 

with  equality,  as  does  the  familiar  coherent  state  |a>  .  Unlike  the  coherent 
state,  which  gives  <Aa^2>  *  <Aa^2>  =  1/4,  (107)  shows  that  there  is  an 

asymmetric  noise  division  between  the  quadratures  (a  noise  squeezing)  in  the 

state  |S;u,v>  ,  with  the  low-noise  quadrature  being  less  noisy  than  a  coherent 
state.  This  noise  reduction  can  be  used,  in  principle,  to  effect  important 
performance  improvements  in  optical  communications  [16]  -  [19]  and  precision 
measurements  [20]  -  [23]. 

As  yet,  there  have  been  no  experimental  observations  of  squeezed  state 
light.  Theoretical  studies,  which  employ  varying  degrees  of  idealization, 
indicate  that  such  states  may  be  generated  by  degenerate  four-wave  mixing  (DFWM) 
[39]  -  [42],  as  well  as  a  number  of  other  nonlinear  optical  processes  [14],  [15], 
[43]  -  [46].  We  are  presently  working  on  a  continuous -wave  DFWM  experiment  using 
homodyne  detection  to  generate  and  verify  the  quadrature  noise  squeezing.  In 
this  experiment,  a  single  frequency-stabilized  laser  will  be  used  to  provide 
all  the  input  beams  to  the  four-wave  mixer,  as  well  as  the  local  oscillators 
for  dual -detector  homodyne  detection.  The  results  of  this  oaper  permit  the 
expected  photocurrent  statistics  for  this  experiment  to  be  derived,  including 
the  effects  of  the  laser's  residual  amplitude  and  phase  fluctuations.  Specifi¬ 
cally,  an  iterated  expectation  approach  is  used,  as  in  Section  III.  The  photo¬ 
current  statistics  are  first  obtained  assuming  the  laser  outout  to  be  a  par¬ 
ticular  coherent  state.  This  entails  a  calculation  of  the  four-wave  mixer 
output  state,  along  the  lines  of  [40],  followed  by  a  calculation  of  the  sort 
oerformed  here  in  Section  1 1 1 A .  To  average  over  the  input  laser  fluctuations, 
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we  assign  to  the  coherent  state  value  for  this  laser  a  classical  probability 
distribution.  We  can  then  proceed  as  in  Section  IIIB,  except  that  the  state 
of  the  signal  field  operator  in  the  homodyne  apparatus  is  now  dependent  on  the 
coherent-state  value  of  the  local -oscillator  field  in  that  apparatus,  because 
both  fields  are  derived  from  the  same  laser. 
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Footnotes 

1.  The  convention  we  use  for  this  Fourier  transform  is  necessitated  by 
the  accepted  quantum-optics  definition  for  what  constitutes  a  positive- 
frequency  field. 

2.  For  a  photomultiplier  tube,  the  internal  current  gain  amplifies  the  current 

we  are  analyzing  by  a  sufficient  amount  to  warrant  its  treatment  as  a  classical 
entity.  In  the  coherent  optical  detection  cases  that  follow,  the  mixing  gain 
produced  by  the  strong  local  oscillator  has  a  similar  effect,  see  [22]. 

3.  A  classical  state  is  either  a  Glauber  coherent  state  or  a  classically  random 
mixture  of  such  states.  In  either  case,  the  density  operator  o  has  a 
proper  P-representation  (17),  The  terminology  arises,  see  below,  because  a 
classical  state  3  gives  rise  to  the  same  statistics  in  quantum  photodetection 
theory  as  found  for  a  classical  field  in  semiclassical  photodetection  theory. 

A.  Because  our  idealized  detector  model  neglects  internal  noise  sources  (dark 
current  shot  noise,  thermal  noise,  etc.)  N  from  Eq.  (21)  corresponds  to  the 
output  of  a  pulse-discriminator/ counter  applied  to  the  output  current 

dx  J(x,t).  In  other  words,  Eq.  (21)  models  the  output  of  an  ideal  (unity 


quantum  efficiency)  photomul tipi ier-tube/pul se-counter  setup. 

Our  choice  for  the  beam  splitter  transformation  agrees  wi-th  that  employed  in 
[18],  and  implies  that  the  field  leaving  the  other  port  of  this  optical  element 
is  -(1-0  E^(x,t)+£  E^g(x,t).  Other  beam  splitter  relations  (see,  e.g. 

[25]  ,  [30])  are  equivalent  to  ours  after  redefinition  of  the  input  and 
output  planes. 

A  critical  aspect  of  the  strong  local  oscillator  condition  acting  through  the 
measurement  operator  (A8)  is  that  the  mean  local  oscillator  field  and  its 
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quantum  noise  both  contribute  to  through  the  direct  detection 

term,  but  only  the  mean  local  oscillator  field  (not  its  quantum  noise) 
contributes  to  JhQ(T)  through  the  mixing  term. 

Very  interesting  noise  reductions  can  accrue  when  the  signal  and  image 
bands  are  quantum-mechanical ly  dependent  [22].. 

8ecause  the  in-phase  and  quadrature  components  of  are  non-commuting 
observables,  the  image  band  noise  enters  into  heterodyning  in  order  to  enforce 
the  Heisenberg  uncertainty  principle  on  ideal  simultaneous  observations  of 
these  incompatible  quantities  (see  [36],  [37]). 

Implicit  in  this  conditioning  statement  is  the  fact  that  the  local  oscillator 
must,  with  very  high  probability,  remain  sufficiently  strong  to  ensure  the 
validity  of  the  Section  II  theory.  Also  note  that  the  signal  field  statistics 
may  depend  on  the  value  of  the  local  oscillator  field,  such  as  occurs  in  a 
laboratory  experiment  when  the  same  laser  is  used  to  obtain  both  the  signal 
and  local  oscillator  beams  (see  Section  IV). 

The  local  oscillator  fluctuations  must  not  be  such  as  to  invalidate  the 
Section  II  theory  for  any  state  i01lq>  that  occurs  with  appreciable  probability. 
Also,  the  signal  state  (density  operator)  may  depend  on  the  value  of  the  local 
oscillator  field,  if,  for  example,  both  beams  originate  from  the  same  laser 
(see  Section  IV) 

For  example,  to  keep  this  added  noise  below  10%  (in  standard  deviation)  of 

-  4 

the  coherent-state  signal  quantum  noise  when  <a^  a^>  =10  ,  we  can  tolerate  no 

more  than  0.3%  local  oscillator  amplitude  fluctuation.  This  limitation  may 
be  significant  in  precision  measurement  applications  for  which  signal-to-noise 
ratios  far  in  excess  of  40dB  are  sought. 
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Fiqure  Captions 

Fig.  1  Geometry  of  an  idealized  surface  photoemitter  of  active  region  A d 
Fig.  2  Configuration  for  optical  homodyne  detection. 

Fig.  3  Configuration  for  optical  heterodyne  detection. 

Fig.  4  Configuration  for  dual-detector  coherent  optical  detection. 
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ABSTRACT 

Degenerate  four-wave  mixing  (DFWM)  lineshapes  are  investigated 
in  sodium  vapor  near  the  0^  resonance  line  using  nearly  Fourier 
transform  limited  pulses.  At  low  pump  intensities  sub-Doppler  reso¬ 
lution  is  obtained.  When  the  Rabi  frequency  associated  with  the 
pump  intensity  becomes  equal  to  the  ground  state  nyperfine  frequency 
separation  of  sodium,  each  component  of  the  double-peaked  DFWM 
spectrum  further  splits  into  two  components  each.  Adiabatic  follow¬ 
ing  model  explains  the  near-resonant  intensity  dependence  of  the 
DFWM  signal  but  is  insufficient  to  explain  the  whole  structure. 
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Resonant  degenerate  four-wave  mixing  (DFWM)  has  been  shown  to  be  an  efficient 
technique  for  the  generation  of  phase  conjugate  wave  fronts^.  In  backward  con¬ 
figurations,  where  the  two  pump  beams  as  well  as  the  two  signal  (probe  and  con¬ 
jugate)  beams  counter-propagate  to  each  other,  the  DFWM  spectrum  is  Doppler  free 
and  has  been  used  for  high  resolution  spectroscopy  .  However,  true  spectral 
information  is  obtained  only  when  the  pump  beam  intensities  are  small  compared 
to  the  saturation  intensity  of  the  resonant  medium.  At  higher  pump  intensities, 
complicated  line  shapes  are  observed;  in  particular,  the  line  broadens  and  splits 
into  two  components.  Several  experimental  and  theoretical  studies  of  this 
behavior  have  appeared  in  recent  literature^.  Most  of  the  studies  have  employed 
cw  lasers  and  hence  steady  state  analyses  have  been  used  to  explain  the  observed 
behavior. 

The  early  experiments  on  resonant  DFWM  employed  pulsed  dye  lasers^.  DFWM 
spectrum  in  this  case  is  also  complicated  and  shows  a  dip  to  zero  at  the  line 
center.  Bloom  et  al .  were  the  first  to  observe  such  a  line  shape.  In  the  pulsed 
case,  fluctuations  (pulse  to  pulse  energy,  pulse-shape,  and  center  frequency 
fluctuations)  and  other  nonlinear  effects  such  as  self  focussing  complicate  the 
line  shape  even  further.  Recently  we  have  measured  the  DFWM  line  shapes  with  a 
pulsed  dye  laser  system  which  was  stabilized  in  center  frequency  and  bandwidth. 

Our  control  of  the  fluctuations  in  the  DFWM  signal  was  essential  for  the  measure¬ 
ment  of  the  photon  counting  statistics  of  light  generated  via  DFWM^.  In  this 
letter  we  describe  the  results  of  a  study  of  DFWM  line  shapes  with  such  a  source. 

The  experimental  setup  is  essentially  the  same  as  that  in  Ref.  6  and  is 
shown  in  Fig.  1.  An  externally  stabilized  cw  ring  dye  laser  (sub  MHz  line  width) 
is  amplified  through  a  chain  of  pulsed  dye  laser  amplifiers  pumped  by  the  smoothed 
output  of  a  Nd:YAG  laser^.  The  output  pulses  of  4  ns  duration  with  typically 
10-20%  energy  fluctuations  and  a  total  line  width  of  200  MHz  which  is  approxi¬ 
mately  twice  the  Fourier  transform  limited  line  width  are  used  to  perform  DFWM. 
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About  10-50  MW/cm  pulse  intensities  are  available  from  such  a  system.  Backward 

Q 

DFWM  geometry  is  employed  with  orthogonally  polarized  pump  beams  .  Phase  conju¬ 
gate  (PC)  signal  whose  polarization  is  orthogonal  to  the  probe  beam  (PB)  is  sepa¬ 
rated  using  a  polarization  beam  splitter  (PBS)  and  directed  onto  a  photomultiplier 
(PMT)  whose  output  is  sent  to  a  boxcar  integrator.  The  output  of  the  boxcar  which 
is  proportional  to  the  PC  pulse  energy  goes  to  the  y-axis  of  a  chart  recorder 
whose  x-axis  is  swept  with  the  dye  laser  frequency.  The  boxcar  whose  gate  dura¬ 
tion  is  chosen  to  be  80  ns  is  triggered  optically.  The  PC  signal  is  first  delayed 
optically  and  then  electrically  using  a  50 n  cable  inserted  between  the  PMT  and 
the  boxcar  in  such  a  way  that  it  arrives  in  the  middle  of  the  gate  duration. 

DFWM  is  performed  in  sodium  vapor  generated  in  a  heat-pipe  oven.  It  is  main¬ 
tained  at  310°C  implying  a  sodium  vapor  density  of  4  x  10^ 4  atoms/cm^  for  the 
measurements  reported  in  this  letter.  2.2  Torr  of  helium  is  used  as  the  buffer 
gas.  All  the  beams  arrive  in  time  coincidence  at  the  sodium  cell.  The  dye  laser 
frequency  is  scanned  over  10  GHz  across  the  sodium  D2  line.  The  spectra  obtained 
in  this  way  are  shown  in  Fig.  2  as  a  function  of  the  pump  beams  intensities  which 
are  kept  the  same  (to  within  10%)  for  both  the  pumps.  Appropriate  neutral  density 
filters  are  introduced  in  the  probe  beam  path  to  avoid  saturation  of  the  PMT.  The 
probe  pulse  energy  is  kept  less  than  a  few  percent  of  the  pump  pulse  energy  in  all 
cases  to  avoid  pump  depletion. 

At  the  lowest  pump  intensities  used  of  0.3  kW/cm2  (all  relative  intensity 

Q 

measurements  are  ±  10%)  ,  the  line  shape  consists  of  two  peaks  which  are  4.5  GHz 
apart  as  shown  in  Fig.  2a.  Also  shown  is  the  fluorescence  observed  in  a  direc¬ 
tion  making  a  small  angle  (=1°)  with  the  PC  beam.  The  pulse  widths  (FWHM)  are 
0.74  and  0.83  GHz  for  the  lower  and  the  higher  frequency  peaks  respectively. 
Therefore  sub-Doppler  resolution  is  obtained  at  these  pump  intensities  even 
though  the  peaks  occur  away  from  the  center  of  the  Doppler-broadened  line  of 


wi dth  1 .8  GHz. 
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As  the  pump  intensities  are  increased  to  0.6  kW/cm  ,  the  higher  frequency  peak 
splits  into  two  components  separated  by  1.3  GHz  as  shown  in  Fig.  2b.  A  further 
increase  of  the  pump  intensities  to  1.9  kW/cm2  leads  to  a  splitting  in  the  lower 
frequency  peak  as  well,  as  shown  in  Figs.  2c  and  2d.  In  Fig.  2d  the  reflectivity 
at  the  highest  frequency  peak  is  0.7X.  The  reason  for  such  a  low  reflectivity  is 
because  a  relatively  large  angle  was  chosen  between  the  pumps  and  the  probe  beam2. 
This  choice  was  dictated  by  the  low  background  requirement  in  our  DFWM  quantum 
noise  measurements  reported  earlier^.  With  smaller  angle  and  higher  pump  intensi¬ 
ties  (at  least  an  order  of  magnitude  higher  than  those  reported  herein)  we  have 
observed  reflectivities  as  large  as  700.  Under  these  conditions  the  conjugate 
pulse  duration  is  significantly  shorter  than  the  probe  pulse  duration.  A  system¬ 
atic  study  of  this  behavior  has  recently  been  reported^. 

2 

At  still  higher  pump  intensities  of  3-10  kW/cm  ,  only  two  peaks  remain  as 
shown  in  Figs.  2e  and  2f.  A  further  increase  of  the  pump  intensities  leads  to  a 
broadening  and  weakening  of  the  lower  frequency  peak  which  is  consistent  with  the 
observations  of  earlier  workers.  A  slight  broadening  and  weakening-  is  already 
observable  in  Fig.  7g.  Jabr  et  al.^  did  not  observe  the  lower  frequency  peak 
with  the  same  choice  of  the  pump  and  probe  beam  polarizations  that  we  employ. 

Their  pump  intensities  were  at  least  an  order  of  magnitude  larger  than  ours.  At 
high  pump  intensities,  r3lf-defocussing  of  the  pump  beams  on  the  lower  frequency 
side  of  the  resonance  causes  a  reduction  of  the  actual  pump  intensities  in  the 
mixing  medium,  which  results  in  a  lowering  of  the  PC  signal. 

The  splitting  of  the  DFWM  li reshape  (Figs.  2b-2d)  was  not  observed  in 
4 

earlier  experiments  .  This  we  believe  is  because  we  have  used  stable  center 
frequency  nearly  Fourier  transform  limited  pulses  for  the  above  measurements^. 
Furthermore,  at  low  intensities  sub-Doppler  lines  whose  widths  are  limited  by 
power  broadening  are  observed. 


The  nonlinearity  responsible  for  DFWM  in  sodium  vapor  is  due  to  resonantly 

4 

enhanced  electronic  Kerr  effect.  Grischkowsky  et  al.  used  the  adiabatic  follow¬ 
ing  (AF)  model  for  a  two  level  atom  under  pulsed  excitation  to  derive  an  expres¬ 
sion  for  the  DFWM  reflectivity  as  a  function  of  the  pump  intensities  and  detuning. 
When  the  AF  conditions  are  satisfied,^  the  third  order  nonlinear  susceptibility 
is  given  by 

(3)  , _ Nu4A2  m 

X  h3(v-v0)3  (1+A2/A2)3/Z 

where  u  is  the  electric  dipole  moment  of  the  two  level  atom,  A  is  the  pump  fields 

amplitudes  (assumed  equal  for  both  the  pumps),  v-vQ  is  the  detuning  of  the  DFWM 

frequency  from  the  atomic  line  center  frequency,  N  is  the  total  number  of  atoms 

2 

per  unit  volume,  and  A$  is  the  normalized  saturation  intensity  given  by 

2  2  2 

A$  =  h|  (\>-vQ)  |/y.  The  DFWM  reflectivity  is  given  by  R  =  tan  <L  =  k  L 
where  k  =  2irvx^3'/2cn,  n  is  the  linear  refractive  index  of  the  medium,  L  is  the 
interaction  length,  and  the  approximation  is  valid  under  conditions  of  weak 
reflectivity. 

As  the  pump  intensities  are  increased,  the  DFWM  signal  saturates  because  of 

the  saturation  of  X^3  .  This  is  verified  for  a  detuning  of  2.3  GHz  where  the  AF 

conditions  are  satisfied  as  shown  in  Fig.  3.  DFWM  reflectivities  were  measured 

from  Fig.  2  at  a  detuning  marked  by  the  arrows.  Dots  are  the  experimental  data 

2 

and  the  solid  line  is  a  fit  to  Eq.  (1)  for  Ig  =  7.8  kW/cm  which  agrees  well  with 

2  2 

I  =  c0c(h(v-v  )/u)‘/2  =  2.1  kW/cm  .  Thus  a  good  agreement  is  obtained  with  the 
AF  moJel  in  its  region  of  validity. 

The  detailed  lineshapes  cannot  be  predicted  using  the  AF  model  as  it  is  not 
valid  sufficiently  close  to  resonance.  Steady  state  models3  developed  to  explain 


cw  DFWM  lineshapes  cannot  be  applied  in  the  transient  case  of  our  experiment. 


The  splitting  of  the  DFWM  lineshape  as  shown  in  Figs.  2b-2d  occurs  when  the  Rabi 
frequency  corresponding  to  the  pump  intensities  is  close  to  the  ground  state 
hyperfine  splitting  of  the  sodium  atom  which  is  1.77  GHz.  Moreover,  the  splitting 
of  the  lower  frequency  peak  occurs  at  approximately  twice  the  pump  intensities 
than  that  of  the  higher  frequency  peak.  This  suggests  that  the  multiple  level 
nature  of  the  sodium  atom  is  playing  a  role.  This  is  not  surprising  because  the 
spectrum  of  the  pulses  used  in  this  experiment  is  much  narrower  than  the  ground 
state  hyperfine  splitting  of  sodium.  For  our  200  MHz,  4  ns  pulses,  sodium  can  be 
well  modelled  as  a  three  level  atom  of  A  type.  The  DFWM  mechanism  in  our  experi¬ 
ment  is  more  complicated  than  it  seems  because  at  low  pump  intensities  of 
Fig.  2a,  although  sub-Doppler  resolution  is  obtained,  the  peaks  are  not  separated 
by  1.77  GHz  as  would  be  the  case  in  a  cw  DFWM  experiment  when  the  spectrum  of  the 
laser  used  is  much  narrower  than  the  ground  state  hyperfine  splitting.  A  detailed 
comparison  with  the  theory  can  only  be  made  when  Maxwell -Bloch  equations  are 
solved  in  the  transient  regime  for  the  backward  DFWM  configuration  taking  the 
three  level  nature  of  the  medium  into  account. 

This  research  was  supported  in  part  by  the  Office  of  Naval  Research. 
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-1 


which  is  satisfied  for  4  ns  pulses  with  =  10  MHz  and  jv-v  |  =2.3  GHz. 
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Figure  Captions: 

1.  Schematic  of  the  experimental  apparatus.  VAG  =  frequency  doubled  NdrYAG 
laser,  DL  =  cw  dye  laser,  M  =  total  reflector,  BS  =  beam  splitter, 

DET  =  detector,  HWP  =  half  wave  plate,  PBS  =  polarization  beam  splitter, 

PMT  =  photo-multiplier  tube,  PC  =  phase  conjugate  beam,  BP  =  backward 
pump  beam,  FP  =  forward  pump  beam,  PB  =  probe  beam. 

2.  DFWM  lineshapes  at  various  pump  intensities.  Vertical  scale  is  arbitrary 

and  is  linearly  proportional  to  the  DFWM  signal,  plotter  scale  factors  are  • 

labelled.  Also  superimposed  are  the  fluorescence  spectra  in  parts  a,  b, 

2 

and  g.  Pump  intensities  are  labelled  in  kW/cm  in  each  plot.  Start  fre¬ 
quency  in  parts  f  and  g  is  slightly  shifted.  9 

3.  Dependence  of  the  DFWM  signal  on  the  pump  intensities. 
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